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ABSTRACT 

Context. Studies of dense molecular-cloud cores at (sub)millimetre wavelengths are needed to understand the early stages of star 
formation. 

Aims. We aim to further constrain the properties and evolutionary stages of dense cores in Orion B9. The prime objective of this study 
is to examine the dust emission of the cores near the peak of their spectral energy distributions, and to determine the degrees of CO 
depletion, deuterium fractionation, and ionisation. 

Methods. The central part of Orion B9 was mapped at 350 pm with APEX/SABOCA. A sample of nine cores in the region were 
observed in C 17 0(2 - 1), H 13 CO + (4 - 3) (towards 3 sources), DCO + (4 - 3), N 2 H + (3 - 2), and N 2 D + (3 - 2) with APEX/SHFI. These 
data are used in conjunction with our previous APEX/LABOCA 870-yum dust continuum data. 

Results. All the LABOCA cores in the region covered by our SABOCA map were detected at 350 pm. The strongest 350 /jm emission 
is seen towards the Class candidate SMM 3. Many of the LABOCA cores show evidence of substructure in the higher-resolution 
SABOCA image. In particular, we report on the discovery of multiple very low-mass condensations in the prestellar core SMM 6. 
Based on the 350-to-870 pm flux density ratios, we determine dust temperatures of r dllst ^ 7.9 - 10.8 K, and dust emissivity indices 
of [3 ~ 0.5 - 1.8. The CO depletion factors are in the range /d ~ 1.6 - 10.8. The degree of deuteration in N 2 H + is ^ 0.04 - 0.99, 
where the highest value (seen towards the prestellar core SMM 1) is, to our knowledge, the most extreme level of N 2 H + deuteration 
reported so far. The level of HCO + deuteration is about 1-2%. The fractional ionisation and cosmic-ray ionisation rate of H2 could be 
determined only towards two sources with the lower limits of ~ 2 - 6 x 10~ 8 and ~ 2.6 x 10~ 17 - 4.8 x 10~ 16 s" 1 , respectively. We also 
detected D2CO towards two sources. 

Conclusions. The detected protostellar cores are classified as Class objects, in agreement with our previous SED results. The 
detection of subcondensations within SMM 6 shows that core fragmentation can already take place during the prestellar phase. The 
origin of this substructure is likely caused by thermal Jeans fragmentation of the elongated parent core. Varying levels of f D and 
deuteration among the cores suggest that they are evolving chemically at different rates. A low fi> value and the presence of gas-phase 
D 2 CO in SMM 1 suggest that the core chemistry is affected by the nearby outflow. The very high N 2 H + deuteration in SMM 1 is 
likely to be remnant of the earlier CO-depleted phase. 

Key words. Astrochemistry - Stars: formation - ISM: abundances - ISM: clouds - ISM: molecules - Radio lines: ISM - Submillimeter: 
ISM 



Using the LABOCA bolometer array on APEX (Atacama 
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New stars in our Galaxy form predominantly in the so-called gi- 
ant molecular clouds (GMCs). The nearest GMC to the Sun is „ , , , r , , , 
, _ . , , , , ' . m _ . , , , Pathfinder Experiment), we mapped that region at 870 urn 

the Orion molecular cloud complejrl This cloud complex con- j A/r . ... . f or >nn i u a n t\ a * ,• 

, . . . r . , f • ■ (Miettine n et al. 20091 hereafter Paper I). The dust continuum 

tains a number of star-forming regions, one example being the — — , . — ~r ,. f. , , r ,. , 

. ., , , J . „ , 1 , f, mapping resulted in the discovery of twelve dense cores of which 

Orion B9 region in the central part of Orion B molecular cloud. c c j . .j-.i_Tr.Ao j 

, ? ? ■ , • , • four were found to be associated with IRAS point sources, and 

5— i Orion B9 represents a very early stage of star formation which is . , . c , , , „i 

rrt ■ -_c ■ . ^ . <; 1 „ ■ ■ , , • eight of them appeared previously unknown sources. Of the 

LU , manifest in the fact that most dense cores in the cloud are ei- , ,. , f , „ , ~ ... 

, . _ „, , . newly discovered cores, two (we called SMM 3 and 4) were 

ther prestellar or contain Class protostars. These objects are „ , . , . , ... „ . _ . , _„ 

r found to be associated with Spitzer 24- and 70-pm sources, 

~Z , a- . . ~ 77 » „■ and were tentatively classified as Class-0 protostellar candidates. 

Send offprint requests to: O. Miettinen J , , . . 

* ti,- ur 1 1 , . . , ... Six cores showed no signs of embedded infrared emission, and 

This publication is based on data acquired with the Atacama b p- , 

Pathfinder Experiment (APEX) under programmes 079.F-9313A, were thus classified as starless. Figure □ shows the 870 ^m 
084.F-9304A, and 084.F-9312A. APEX is a collaboration between emission as contours overlaid on the Spitzer/MlPS 24-p.m lm- 
the Max-Planck-Institut fur Radioastronomie, the European Southern age of the central part of Orion B93. In Paper I, we also de- 
Observatory, and the Onsala Space Observatory. 

** Our SABOCA map shown in Fig. 2 is available electronically in 2 Comparison of our LABOCA map to the new SABOCA map to- 

FITS format at the CDS. gether with Spitzer images showed that the LABOCA observations 

1 In this paper we adopt a distance of 450 pc to the Orion GMC were not well-pointed. The SABOCA peaks appeared to systemati- 

(Genzel & Stutzki 1989). The actual distance may be somewhat smaller cally lie southeast from the LABOCA peak positions. The difference 

as, for example, Menten et al. (2007) determined the distance to the in angular resolution between SABOCA and LABOCA data did not 

Orion Nebula to be 414 ± 7 pc. appear to be the cause of the offset in peak positions; by smoothing 
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rived the degree of deuteration, or the /V(N 2 D + )//V'(N 2 H + ) co- 
lumn density ratio, towards selected positions and found the 
values in the range 0.03-0.04, comparable to those seen in 



other low-mass star-forming regions (e.g., Crapsi et al. 2005 
Emprechtinger et al. 2009| IFriesen et al. 20 1 Obl l . Taking advan 



tage of the H 2 D + data from Harju et al. (2006), the ionisation 
degree and the cosmic-ray ionisation rate of H 2 towards the 
same target positions were estimated to be x(e) ~ 10~ 7 and 
f H , ~ 1 - 2 x 10~ 16 s" 1 , respectively. 

The physical properties of dense cores in Orion B9 were 
studied further by Miettinen et al. (2010; hereafter Paper II). 
Using the observations of the NH3 (1,1) and (2, 2) inversion 
transitions performed with the Effelsberg 100-m telescope, we 
determined the gas kinetic temperature, kinematical properties, 
and dynamical state of the cores. The gas kinetic temperature of 
the cores were found to be 7\j n ~ 9.4 - 13.9 K, and the inter- 
nal non-thermal motions in the cores appeared to be subsonic, or 
at most transonic. The virial-parameter analysis showed that the 
starless cores in the region are likely to be gravitationally bound, 
and thus prestellar objects. Interestingly, some of the cores were 
found to have significantly (by ~ 5 — 7.5 km s~') lower ra- 
dial velocity than the systemic velocity of the region (vlsr ~ 9 
km s" 1 ). This suggests that they belong to the "low-velocity part" 
of Orion B which is likely to originate from the feedback from 
the massive stars of the nearby Ori OB lb association. 

In this paper, we present the results of our new APEX 
observations of the Orion B9 region. These include the dust 
continuum mapping with the SABOCA (Submillimetre APEX 
Bolometer Camera) bolometer array at 350 /mi, and observations 
of the molecular-line transitions C 17 0(2 - 1), H 13 CO + (4 - 3), 
DCO + (4 - 3), N 2 H + (3 - 2), and N 2 D + (3 - 2). The SABOCA 
350-/im observations allow us to probe the peak of the core's 
spectral energy distributions (SEDs). We also aim to study the 
degrees of CO depletion and deuterium fractionation, and the 
fractional ionisation in the cores. These are useful parameters 
to further constrain the physical and chemical properties of the 
sources, and their evolutionary stages. 

This paper is structured as follows. The observations and 
data-reduction procedures are described in Sect. 2. The direct 
observational results are presented in Sect. 3. In Sect. 4, we des- 
cribe the analysis and present the results of the physical and 
chemical properties of the cores. Discussion of our results is pre- 
sented in Sect. 5, and in Sect. 6, we summarise the main conclu- 
sions of this study. 



2. Observations and data reduction 

2.1. Submillimetre dust continuum observations 

The central 1 8 f0 X 14!5 part of Orion B9 was mapped with 
SABOCA (|Siringo et al. 2010} on the APEX 12-m telescope 
(Giis ten et al. 20061 1 at Llano de Chajnantor (Chilean Andes). 
SABOCA is a 37-channel on-sky bolometer array operating at 
350 /mi, with a nominal resolution of ~ 7'.'5 (HPBW). The ef- 
fective field of view of the array is 1 '5. The SABOCA passband 



the SABOCA map to correspond the resolution of our LABOCA data, 
and regridding the maps onto the same grid, the offsets still remained. 
We used the brightest SABOCA 350-/jm source in the region, SMM3, 
as a reference to adjust the pointing of our LABOCA map. The point- 
ing of the LABOCA map was shifted by (Aa, AS) = (+1'.'72, -9T97). 
Unfortunately, the target positions of our molecular-line observations, 
which were chosen to be the LABOCA peak positions from our previ- 
ous map, are now slightly offset from the dust maxima (but still within 
the beam size for most line observations). 



has an equivalent width of about 120 GHz centred on an effective 
frequency of 852 GHz. 

The observations were carried out on 5 October and 10 
November 2009. The atmospheric zenith opacity at 350 /mi was 
measured using the sky-dip method, and was found to be in 
the range Tz 50/Jm = 0.8 - 1.1. The amount of precipitable wa- 
ter vapour (PWV) was in the range 0.4-0.6 mm. The telescope 
pointing and focus checks were made at regular intervals using 
the planets Mars, Jupiter, and Uranus, and several secondary ca- 
librators (e.g., CRL618, N2071IR, and CW Leo). The absolute 
calibration uncertainty for SABOCA is 25-30%. Mapping was 
performed using a 'fast-scanning' method without chopping the 
secondary mirror (Reic hertz et al. 20011 >. The mosaic was con- 
structed by combining 147 individual fast-scanning (typically 
Y.I s~') maps. The total integration time spent on the area was 
9.7 hr. 

Data reduction was done with the CRUSH-2 
(Comprehensive Reduction Utility for SHARC-2) (version 
2.03-2) software package (Kovacs 2008), which includes cali- 
bration data covering the period of our SABOCA observations. 
We used the pipeline iterations with specifying the 'deep' op- 
tion, which is appropriate for point-like sources. To improve the 
image appearance and increase the signal-to-noise (S/N) ratio, 
a beam-smoothing was applied, i.e., the map was smoothed 
with a Gaussian kernel of the size 7'.'5 (FWHM). Therefore, the 
angular resolution of the final image is 1076 (0.02 pc at 450 pc). 
The gridding was done with a cell size of 1'.'5. The resulting lcr 
rms noise level in the final co-added map is ~ 0.06 Jy beam 1 . 



2.2. Molecular-line observations 

The spectral-line observations of C n O(2 - 1), DCO + (4 - 3), 
N 2 H + (3 - 2), and N 2 D + (3 - 2) towards nine Orion B9 cores, 
and H 13 CO + (4 - 3) observations towards three cores, were car- 
ried out on 2-3, 8-9, and 11-12 September, and 13 and 18-19 
November 2009 with APEX. The target cores with their physi- 
cal properties derived in Paper II are listed in Table Q] In Table Q] 
we give the LABOCA submm peak positions of the cores after 
adjusting the pointing. The target positions of our molecular-line 
observations can be found in Table 1 of Paper II. The observed 
transitions, their spectroscopic properties, and observational pa- 
rameters are listed in Tablef2] The N 2 H + (3-2) observations were 
already presented in Paper II, but are re-analysed in the present 
paper. 

As frontend for the C n O(2- 1) and N 2 D + (3-2) observations, 
we used APEX-1 of the SHFI (Swedish Heterodyne Facility 

IVassilev etal. 2008al b). 



Instrument; Belitsky et al. 2007 



APEX-1 operates in single-sideband (SSB) mode using side- 
band separation mixers, and it has a sideband rejection ratio 
> 10 dB. The SHFI sidebands are separated by 12 GHz, i.e., 
±6 GHz around the local oscillator (LO) frequency. Therefore, 
the centre frequency for the image band is 12 GHz above 
or below the observing frequency depending on whether the 
receiver tuning is optimised for operation in the upper or lower 
sideband (USB or LSB), respectively. For the H 13 CO + (4 - 3), 
DCO + (4 - 3), and N 2 H + (3 - 2) observations the frontend used 
was APEX-2, which has similar characteristics as APEX-1. The 
backend for all observations was the Fast Fourier Transfrom 
Spectrometer (FFTS; IKlein et al.~2 006) with two 1 GHz units. 
Both units were divided into 8 192 channels. The two units 
were connected to one receiver each, thus providing 1 GHz 
bandwidth for two receivers simultaneously. 
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Fig. 1. Spitzer/MTPS 24-ji/m image of the central part of Orion B9 overlaid with contours showing the LABOCA 870-/im dust 
continuum emission. The contours go from 0.1 (~ 3.3a) to 1.0 Jy beam -1 in steps of 0.1 Jy beam -1 . The 24-fim image is shown 
with a logarithmic scaling to improve the contrast between bright and faint features. The colour bar indicates the 24-yum intensity 
scale in units of MJy sr -1 . The green plus signs show the target positions of our molecular-line observations (i.e., the submm peak 
positions of the LABOCA map before adjusting the pointing; see text). The green line shows the base and northwest tip of the HH 
92 jet driven by IRAS 05399-0121 ( |Bally et al. 2002} . The 0.2-pc scale bar and the effective LABOCA beam HPBW (~ 20") are 
shown in the bottom left (Paper I). 



The observations were performed in the wobbler-switching 
mode with a 100" azimuthal throw (symmetric offsets) and a 
chopping rate of 0.5 Hz. The telescope pointing and focus cor- 
rections were checked by continuum scans on the planets Mars 
and Uranus, and the pointing was found to be accurate to ~ 3". 
Calibration was done by the chopper-wheel technique, and the 
output intensity scale given by the system is 7^, the antenna 
temperature corrected for atmospheric attenuation. The observed 
intensities were converted to the main-beam brightness tem- 
perature scale by r M B = T^/t/mb, where 77MB = B e g/F eS is the 
main beam efficiency, and B e s and F e ff are the beam and forward 
efficiencies, respectively. The absolute calibration uncertainty is 
estimated to be 10%. 

The spectra were reduced using the CLASS90 programme of 
the IRAM's GILDAS software packag^- The individual spectra 
were averaged and the resulting spectra were Hanning-smoothed 
in order to improve the S/N ratio of the data. A first- or third- 
order polynomial was applied to correct the baseline in the final 
spectra. The resulting Icr rms noise levels are ~ 20 - 80 mK at 
the smoothed resolutions. As shown in the second last column of 
Table [2] the on-source integration times were different between 
different sources. This is particularly the case for C 17 0(2 - 1), 



N 2 D + (3 - 2), and H 13 CO + (4 - 3) where f int varies by a factor of 
~ 3 — 7. This explains the variations in the rms noise values (up 
to a factor of 2.4 in the C n O data). 

The 7 = 2-1 transition of C 17 contains nine hyperfine 
(hf) components. We fitted this hf structure using "method hfs" 
of CLASS90 to derive the LSR velocity (vlsr) of the emission, 
and FWHM linewidth (Av). The hf line fitting can also be used 
to derive the line optical thickness, r. However, in all spectra the 
hf components are mostly blended together and thus the opti- 
cal thickness could not be reliably determined. For the rest fre- 
quencies of the hf components, we used the values from Ladd 
et al. (1998; Table 6 therein). The adopted central frequency, 
224714.199 MHz, is that of the J F = 2 9/2 -> I7/2 hf compo- 
nent which has a relative intensity of = i. 

Also, rotational lines of H 1 3 CO + and DCO + also 
have hf structure (e.g., Schmid-Burgk et al. 2004 



http : //www . iram . fr/IRAMFR/GILDAS 



ICaselli & Pore 20051 1. The J = 4 - 3 transition of DCO + is split 
up into six hf components. To fit this hf structure, we used the 
rest frequencies from the CDMS database (Mii ller et al. 20051 . 
The adopted central frequency of DCO + (4 - 3), 288 143.855 
MHz, is that of the Jf = 4s — > 34 hf component which has 
a relative intensity of Ri = We note that this frequency 
is 2.7 kHz lower than the value determined by Caselli & 
Dore (2005; their Table 5). The frequency interval of the hf 
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components for DCO + (4 - 3) is very small. Therefore, the 
lines overlap significantly which causes the hf structure to be 
heavily blended. In general, for the 7 U p per > 3 lines of DCO + 
the hf components are so heavily blended that, even determined 
through a single Gaussian fit, the linewidth is not expected to 
be significantly overestimated (Casel fi~& Pore 20051 1. The hf 
structure of H I3 CO + lines is more complicated because both 
the 'H and 13 C nuclei have a nuclear spin of / = 1/2, and the 
nuclear magnetic spin couples to rotation. To our knowledge, 
the rest frequencies of the H 13 CO + (4 - 3) hf components 
have not been published. We thus fitted these lines using a 
single Gaussian fit to derive the values of Vlsr and Av. As in 
the case of DCO + (y up pei > 3), the Av thus determined is not 
expected to be significantly overestimated because of the strong 
blending of the hf components. The central frequency used 
was 346998.344 MHz (CDMS), which is 3.0 kHz lower than 
the value determined by Schmid-Burgk et al. (2004; Table 3 
therein). 

The 7 = 3-2 transitions of both N 2 H + and N 2 D + contain 38 
hf components. The hf lines were fitted using the rest frequencies 
from Pagani et al. (2009b; Tables 4 and 10 therein). The adopted 
central frequencies of N 2 H + (3 - 2) and N 2 D + (3 -2), 279 5 1 1 .832 
and 231 321.912 MHz, are those of the J FlF = 3 45 -> 2 34 hf 
component which has a relative intensity of R, = ||. Also in 
these cases, the hf components are blended and thus the value of 
t could not be reliably determined through hf fitting. 



3. Observational results 

3.1. SABOCA 350-yim emission 

The 350-/im SABOCA map is shown in Fig. |2] Almost all the 
cores detected with LABOCA show also clear 350 fim emission. 
The exceptions are SMM 5 and Ori B9 N for which the 350-^/m 
peak flux densities are at the levels 3.5<x (210 mJy) and 3<x (180 
mJy), respectively. With the peak flux density of 60.5<x (3630 
mJy), SMM 3 is by far the strongest 350-^/m source in the re- 
gion. Four of the cores are resolved into at least two emission 
peaks in the SABOCA image. In the case of SMM 3, there is a 
350-^m condensation, we call SMM 3b (4.2<x or 250 mJy), at 
about 36" east of the "main" source. At about 17" from SMM 
3b, there is another Act (240 mJy) emission peak, designated 
here as SMM 3c. The 870 /mi emission of SMM 3 extends to 
the direction of the subcondensations SMM 3b and 3c, and both 
of them lie either within or at the borderline of the 3.3<x (99 
mJy) 870-/im contour. SMM 4 is resolved into two condensa- 
tions. The western condensation coincides with the LABOCA 
peak position, whereas the eastern one, SMM 4b, is coincident 
with a Spitzer 24-/mi source near SMM 4 (see Sect. 5.7.2). 
The elongated core SMM 6 is resolved into at least three sub- 
condensations. The northwesternmost condensation corresponds 
to LABOCA peak of SMM 6. Finally, SMM 7 shows a hint 
of substructure in the western part of the core, where the 870 
//m emission is above 3.3<r (99 mJy) level. SMM 3, 4, 6, 7, 
and their substructure will be discussed further in Sect. 5.7. 
There are also two 4 - 4.5<x (240 - 270 mJy) emission peaks at 
»20oo.o = 05 h 42 m 47.4 s ,5 20 oo.o = -01°17'15" (south of SMM 3b) 
anda 20 oo.o = 05 h 43 m 24.7 s , 5 200 o.o = -01°14'35" (southeast of 
SMM 7), and a few 3.5 -4cr (210 - 240 mJy) peaks to the east of 
SMM 6 at a 200 o.o = 05 h 43 m 16.5 s , <5 20 oo.o - -01°19'00". There 
is some weak LABOCA emission just slightly north of the latter 
peaks. However, because these "additional" sources are below 
5cr (300 mJy), were not detected by LABOCA at 870 /mi, and 



were not identified by the SIMBAD Astronomical Databas^J 
they may well be unreal and are not discussed further in this 
paper. 

To identify and extract the cores from the SABOCA map, 
we employed the commonly used two-dimensional clumpf ind 
algorithm, clfind2d, developed by Williams et al. (1994). The 
clfind2d routine determines the peak position, the FWHM size 
(not corrected for beam size), and the peak and total integrated 
flux density of the source based on specified contour levels. The 
algorithm requires two configuration parameters: i) the intensity 
threshold, i.e., the lowest contour level, T\ ow , which determines 
the minimum emission to be included into the source; and it) the 
contour level spacing, AT, which determines the required "con- 
trast" between two sources to be considered as different objects. 
We set both parameters to 3<x (180 mJy). The selected 3cr con- 
tour levels turned out to give the best agreement with the identi- 
fication by eye. With these parameter settings, the SMM 6c and 
6d condensations are treated as a single source by clumpfind. 

The J2000.0 coordinates of the peak 350 yum emission, 
source effective radius (R e s = ~\/A/n, where A is the projected 
area within the 3<x contour), and peak and integrated flux den- 
sities are listed in Cols. (2)-(6) of Table [3] In Col. (7), we also 
list the flux densities measured in a 40" diameter aperture from 
the SABOCA map smoothed to the resolution of our LABOCA 
data. These flux densities were determined using the CRUSH-2 
programme. The effective radius is only given for sources which 
are larger than the beam size. The total flux density uncertainty 

was derived from <t(Sa) = + °"s> where cr ca i is the abso- 

lute calibration error (adopted to be 30% of flux density), and 
<xs is the uncertainty in the flux-density determination based on 
the rms noise near the source area. The lcr rms noise in the 
smoothed SABOCA map, <x = 0.08 Jy beam 1 , is slightly higher 
than in the original map. We note that the data reduction pro- 
duces negative artefacts ("holes") around regions of bright emis- 
sion, most notably around SMM 3. This decreases the source's 
peak intensity, introducing an additional uncertainty in the flux 
density. The uncertainties due to negative bowls are neglected in 
the subsequent analysis. 

3.2. Spectra 

The Hanning-smoothed spectra are shown in Fig. [3] The 
C 17 0(2 - 1) line is clearly detected towards all sources except 
SMM 4 where only a low-velocity component is observed at 
vlsr - 1.7 km s _1 (hereafter, SMM 4-LVC, where LVC stands 
for low-velocity component). We note that SMM 4 shows strong 
NH 3 (1, 1) emission at about 9.1 km s" 1 (Paper II). IRAS 05405- 
01 17 (henceforth, IRAS05405 etc.) shows two additional veloc- 
ity components at ~ 1.3 and 3.0 km s , and Ori B9 N shows an 
additional and wide C n O(2 - 1) line at ~ 1.9 km s _1 . There is 
also a hint of C 17 0(2- 1) emission at the systemic velocity of the 
region (~ 9 km s _1 ) in the spectrum towards SMM 7, although 
the radial velocity of SMM 7, as determined from NH3(1, 1) 
measurements in Paper II, is 3.6 km s . The hf-structure of the 
C n O(2 - 1) line is partially resolved in IRAS05399, SMM 1 and 
3, and Ori B9 N. 

The H 13 CO + (4 - 3) observations were carried out only to- 
wards three sources (IRAS05399, SMM 1, and SMM 4). The 
line is only detected in IRAS05399 and SMM 4-LVC. The 
DCO + (4 - 3) line is detected towards all the other sources ex- 
cept IRAS05405, SMM 4, and SMM 5. The line is also quite 
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Table 1. Source list. 



Source 


c^ooo.o" 


^2000.0° 




M 




Af(H 2 )* 




<n(H 2 )> 


Class 




[h:m:s] 


[":':"] 


[K] 


[Mo] 




[10 22 cm 


2 ] 


[10 4 cnr 3 ] 




IRAS 05399-0121 


05 42 27.5 


-01 20 00 


13.5 ± 1.6 


7.8 ±0.5 


4.2 


± 0.9/2.8 


± 0.6 


5.5 ± 1.3 


0/1 


C A H A It 1 

SMM 1 


05 42 30.5 


-01 20 55 


11.9 ± 0.9 


11.1 ± 0.3 


2.8 


± 0.4/2.7 


± 0.4 


5.3 ± 0.9 


prestellar 


SMM3 


05 42 45.2 


-01 16 13 


11.3 ± 0.8 


7.8 ± 0.6 


8.4 


± 1.1/2.5 


±0.3 


10.5 ±2.1 





IRAS 05405-0117 


05 43 02.7 


-01 16 31 


11.3 ± 0.6 


2.8 ± 0.4 


1.4 


±0.1/1.2 


±0.1 


3.8 ±0.5 





SMM4 


05 43 04.0 


-01 15 54 


13.9 ±0.8 


2.8 ±0.3 


1.4 


±0.1/1.1 


±0.1 


3.8 ±0.4 





SMM5 


05 43 04.6 


-01 17 17 


11.3 ±0.7 


1.9 ±0.4 


1.2 


±0.1/0.9 


±0.1 


2.5 ± 0.5 


prestellar ? 


SMM6 


05 43 05.2 


-01 18 48 


11.0 ±0.4 


8.2 ± 1.1 


2.0 


±0.1/1.7 


±0.1 


2.2 ±0.3 


prestellar 


Ori B9 N 


05 43 05.7 


-01 14 51 


13.4 ± 1.3 


2.3 ±0.4 


0.9 


±0.1/0.9 


±0.1 


2.1 ±0.4 


prestellar ? 


SMM7 


05 43 22.2 


-01 13 56 


9.4 ± 1.1 


3.6 ± 1.0 


3.4 


±0.9/2.1 


±0.5 


4.8 ± 1.3 


prestellar 



Notes. Columns (2) and (3) give the equatorial coordinates [(», <5) 2 ooo.o]- Columns (4)-(7) list the gas kinetic temperature, core mass, beam- 
averaged peak H 2 column density, and the volume-averaged H 2 number density, respectively. In the last column we give the comments on the 
source classification. The virial parameter of the starless cores SMM 5 and Ori B9 N is a v j r ^ 2, and thus it is unclear if they are prestellar (Paper 
II). <D) These coordinates refer to the LABOCA peak positions after adjusting the pointing by (Ace, AS) = (+1'.'72, —9! 91). The coordinates of 
our molecular-line observation target positions can be found in Table 1 of Paper II. <b> The first value is calculated towards the revised LABOCA 
peak position by using the temperature derived towards the line observation target position. The latter N(H 2 ) value refers to the position used 
in molecular-line observations. 



Table 2. Observed spectral-line transitions and observational parameters. 



Transition 


V 


E u /k B 


^crit 


HPBW 


'/MB 


T 

± sys 


PWV 


Channel 


spacing" 




rms 




[MHz] 


[K] 


[cm" 3 ] 


["] 




[K] 


[mm] 


[kHz] 


[km s~'] 


[min] 


[mK] 


CO(2- 1) 


224714.199" 


16.2 


9.0 x 10 3 


27.8 


0.75 


255-288 


0.8-1.4 


122.07 


0.16 


2.9-19.5 


24-58 


N 2 D + (3 - 2) 


231 321.912 r 


22.2 


1.7 x 10 6 


27.0 


0.75 


235-335 


0.2-3.8 


122.07 


0.16 


11.5-37.5 


22-35 


N 2 H + (3 - 2) 


279 511.832 c 


26.8 


2.9 x 10 6 


22.3 


0.74 


211-359 


0.04-1.4 


122.07 


0.13 


2.7-5.5 


32-84 


DCO + (4 - 3) 


288 143.855'' 


34.6 


4.7 x 10 6 


21.7 


0.74 


170-172 


0.01-0.5 


122.07 


0.13 


4.1-6.5 


23-29 


H 13 CO + (4-3) 


346 998.344 f 


41.6 


8.2 x 10 6 


18.0 


0.73 


362-521 


0.9-1.3 


122.07 


0.11 


8.5-22 


25-55 



Notes. Columns (2)-(4) give the rest frequencies of the observed transitions (v), their upper-state energies (E u /k B , where k B is the Boltzmann 
constant), and critical densities. Critical densities were calculated at T = 10 K using the collisional-rate data available in the Leiden Atomic 
and Molecular Database (LAMDA; http : //www. strw. leidenuniv.nl/~moldata/ I (Schoier et al. 2005 1. For N 2 D + , we used the Einstein 
A-coefficient from Pagani et al. (2009b) and the same collisional rate as for N 2 H + . Columns (5)— (12) give the APEX beamsize (HPBW) and 
the main beam efficiency (77MB ) at the observed frequencies, the SSB system temperatures during the observations (r sys in Tmb scale, see text), 
the amount of PWV, channel widths (both in kHz and km s~') of the original data, the on-source integration times per position (/;„,), and the 
ltr rms noise at the smoothed resolution. 

^ The original channel spacings. The final spectra were Hanning-smoothed which divides the number of chan- 
nels by two. (w From Ladd et al. (1998). {c) From Pagani et al. (2009b). (</) From the Jet Propulsion Laboratory (JPL; 
http://spec.jpl.nasa.gov/) spectroscopic database (Pickett et al. 1998). w From the Cologne Database for Molecular Spectroscopy 
(CDMS; http : //www . astro . uni -koeln . de/cdms/catalog) ( Miil ler et al. 2005b . 



weak (Act) in Ori B9 N. Again, clear emission from an addi- 
tional velocity component can be seen towards SMM 4 and Ori 
B9N. 

As already presented in Paper II, N2H + (3 - 2) emission is 
clearly seen in all sources with the exception of target position 
near Ori B9 N where the line is very weak. It should be noted, 
however, that towards SMM 4 only the LVC is detected. Ori B9 
N shows an additional velocity component at about 1.8 km s _1 . 
Finally, the N 2 D + (3 - 2) line was detected in IRAS05399, SMM 
1, SMM 3, IRAS05405, SMM 4-LVC, and SMM 6. In addition, 
the line appears at ~ 3 km s towards IRAS05405. 



3.2.1 . Other line detections 

All the observed sources except SMM 5 show additional spectral 
lines in the frequency band covering the N2D + (3 - 2) transition. 
The line identification was done by using Weeds, which is an 
extension of CLASS (IMaret et al. 201 11) , and the JPL and CDMS 
spectroscopic databases. We used the LTE modelling application 
of Weeds to check if all predicted lines of a candidate molecule 
are present in the observed spectrum. In some cases, we were 



able to reject some line candidates on the basis of non-detection 
of other transitions expected at nearby frequencies. 

The Jk„.k c = 4o,4 - 3o,3 transition of ortho-^CO at ~ 231.4 
GHz was detected towards IRAS05399 and SMM 1 in the LSB 
(see Fig. |H). Also, the J = 2 - 1 transition of C ls O at ~ 219.56 
GHz (LSB) was detected in the image sideband towards all 
sources except SMM 5 [marked with "(i)"; Fig. 0. An addi- 
tional velocity component is detected towards IRAS05405 and 
SMM 4. Moreover, in the spectrum towards Ori B9 N several 
distinct velocity components are detected. The possible fourth 
velocity component of C ls O could, in principle, be blended 
with the line at ~ 231.06 GHz, which can be assigned to ei- 
ther OCS(7 = 19 - 18) (231060.9830 MHz, E u /k B = 110.9 
K) or CH 3 NH 2 -E(./^ = 7 2 , 5 - 7 1>5 ) (231060.6041 MHz, 
Eu/kB = 75.6 K) because these two transitions are also blended. 
However, these can almost certainly be excluded because of the 
high transition energies involved, and because both species are 
expected to be formed by grain-surface chemistry. Detection of 
these species in a cold core is therefore very unlikely. Note that 
the lines "leaking" from the rejected image band are heavily at- 
tenuated by the sideband filter. Therefore, we cannot establish 
the correct intensity scale for the C 18 0(2 - 1) lines. 
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.2 

C 

1 
o 



-0.33 



10.0 05,0 5:43;fl0.0 

Right ascension ( J2000) 




■0.29 



-0.17 



0.03 



0.31 



0.66 



1.1 



1.6 



2.2 



2.9 



3.6 



Fig. 2. SABOCA 350-/mi image (smoothed to 1(X'6 resolution) of the central part of Orion B9 (colour scale and green contours). 
The image is shown with a square root scaling, and the colour bar indicates the flux density in units of Jy beam 1 . The rms level is 
0.06 Jy beam" 1 (1(f). The first 350-/im contour and the separation between contours is 3cr. The white contours show the LABOCA 
870-pm dust continuum emission as in Fig.Q] The yellow plus signs show the target positions of our molecular-line observations. 
The 0.2-pc scale bar and beam HPBW are shown in the bottom right. 



Table 3. The 350-/vm core properties. 





Peak position R eS 


~ peak 
350 


S 350 


o40 a 
350 


Name 


»2ooo.o [h:m:s] fcooo.o [":':"] ["] 


[Jy beam ] 


[Jy] 


[Jy] 



SMM3 


05 


42 


45 


3 


-01 


16 


16 


13.3 


3 


63 


5.4 


± 1.6 


3 


9i 


= 1 


2 


SMM 3b 


05 


42 


47 


6 


-01 


16 


24 







25 














SMM3c 


05 


42 


48 


6 


-01 


16 


32 







24 














IRAS 05405-0117 


05 


43 


03 





-01 


16 


31 


7.6 





47 


0.5 


±0.2 





6 = 


= 


2 


SMM 4 


05 


43 


04 





-01 


15 


52 


5.8 





25 


0.2 


±0.1 





li 


= 


3 


SMM 4b 


05 


43 


05 


7 


-01 


15 


57 


7.5 





49 


0.5 


±0.2 










SMM 5 


05 


43 


03 


8 


-01 


16 


59 







21 









5 = 


= 


2 


SMM 6 


05 


43 


04 


9 


-01 


18 


49 


8.2 





33 


0.5 


±0.2 





7 = 


= 


3 


SMM 6b 


05 


43 


06 


7 


-01 


19 


02 


7.0 





27 


0.3 


±0.1 










SMM 6c 


05 


43 


07 


9 


-01 


19 


13 


6.2* 





21 


0.2: 


tO.l 6 










SMM 6d 


05 


43 


08 


6 


-01 


19 


22 







23 














Ori B9 N 


05 


43 


05 


6 


-01 


14 


55 







18 









5i 


= 


2 


SMM 7 


05 


43 


22 


7 


-01 


13 


55 


10.1 





32 


0.7 


±0.2 





8: 


= 


3 


SMM 7b 


05 


43 


21 


3 


-01 


14 


02 







21 















Notes. (fl) Flux density measured in a 40" diameter aperture from the SABOCA map smoothed to the resolution of the LABOCA map. (i) These 
values refer to the combined size and flux density of SMM 6c and 6d. 



In Table [4] we list the detected extra transitions, their rest 
frequencies, and upper-state energies. The rest frequencies and 
upper-state energies were assigned using JPL, CDMS, and 
Splataloguj^(Remij an et al. 2007[ ) spectroscopic databases. The 
D2CO detections are discussed further in Sect. 5.6. 



http : //www . splatalogue . net/ 



3.3. Spectral-line parameters 

The spectral line parameters are given in Table|5] In this table we 
give the radial velocity (vlsr), FWHM linewidth (Av), peak in- 
tensity (Tub), integrated intensity ( J TmbcIv), peak optical thick- 
ness of the line (to), and excitation temperature (r ex ). For non- 
detections, the 3<t upper limit on the line intensity is given. The 
values of v LSR and Av for C 17 0, DCO + , N 2 H + , and N 2 D + were 
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IRAS 05399-012 
N 2 H + (3-2) 



IRAS 05399-0121 
DCO + (4-3) 

_ 


1 








IRAS 05399-0121 J 
N 2 D + (3-2) 









SMM 1 ft 
DCO + (4-3) 

r\ . .* , i,. a A A A a /-A ) 






5 '0 15 2 
%a [km s "'] 


SMM 1 i 
N 2 H + (3-2) 








SMM 1 

N 2 D + (3-2) | 










SMM 3 
N 2 D + (3-2) 











Fig. 3. Smoothed C 17 0(2 - 1), H 13 CO + (4 - 3), DCO + (4 - 3), N 2 H + (3 - 2), and N 2 D + (3 - 2) spectra. Overlaid on the C n O, DCO + , 
N 2 H + , and N 2 D + spectra are the hf-structure fits. The relative velocities of individual hf components in these transitions are labelled 
with a short bar on the spectra towards IRAS05399. The H 13 CO + spectra are overlaid with a single Gaussian fit. The fits to the lines 
at the systemic velocity ~ 9 km s are shown as green lines, whereas the red lines show fits to the lines at lower velocities. Note 
that H 13 CO + (4 - 3) observations were carried out only towards three sources and the line is detected only towards IRAS05399 and 
SMM 4 (at- 1.5 km s- 1 ). 



Table 4. Other detected species/transitions. 



Species/ 




V 


E u /k B 


transition" 




[MHz] 


[K] 


C ls O(7 = 2- 


If 


219 560.3600 (JPL) 


15.8 


ortho-D 2 CO(J Ka ,K c = 


4o,4 - 3o,3) 


231410.234 (CDMS) 


27.9 



Notes. (a) For asymmetric top molecules such as D2CO, K a and K L 
refer to the projection of the angular momentum along the a and c 
principal axes. The detected 0-D2CO transition exbits a-type selection 
rules (AK a = 0, ±2, . . . and AK C = ±1, ±3, . . .) ( |Gordy & Cook 1984) . 
<b) Seen in the image band. 



derived through fitting the hf structure. For the other lines these 
parameters were derived by fitting a single Gaussian to the line 
profile. Also, the values of Tme and J TMBdv were determined 
from Gaussian fits. The integrated C I7 0(2 - 1) line intensities 
of IRAS05399, SMM 1, 3, and 7, and Ori B9 N were calculated 
over the velocity range given in Col. (6) of Table |5]because the 
hf structure of the line is partially resolved in these cases. The 
uncertainties reported in Vlsr and Av represent the formal lcr 
errors determined by the fitting routine, whereas those in Tmb 
and Jr M Bdv also include the 10% calibration uncertainty. We 
used RADEjfl ( van der Tak et al. 20071 1 to determine the values of 
T and T ex for the lines of C n O, H 13 CO + , DCO + , and N 2 H + . 



6 http : //www . strw . leidenuniv . nl/~moldata/ radex . html 
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RADEX modelling is described in more detail in Sect. 4.4 where 
we discuss the determination of molecular column densities. For 
the rest of the lines, to was determined through Weeds modelling 
(Sect. 4.4), and the associated error was estimated from the 10% 
calibration uncertainty. The lines appear to be optically thin in 
most cases, except that the N2H + lines are mostly optically thick 
(to > 1). The r ex values for C n O(2 - 1) are close to 7\i n , indicat- 
ing that the lines are nearly thermalised. The r ex [H 13 CO + (4-3)] 
values are found to be similar to r ex [DCO + (4 - 3)]. We note that 
the r ex values obtained for the 7 = 3-2 transition of N2FF 
(3.8-5.0 K) are mostly comparable to the values T ex = 4.6 - 5.5 
K determined or assumed in Paper II. For N2D + (3 - 2) we adopt 
as T ex the value derived for N2H + (3 - 2). 



4. Analysis and results 

4. 1 . Core properties derived from 350 fim emission 

We used the 350-yum dust continuum data to determine the 
mass, peak beam-averaged column density of H2, and volume- 
averaged H2 number density of the cores and their condensa- 
tions. The formulas for the mass and H2 column density can be 
found in Paper I [Eqs. (2) and (3) therein]. Note that in Eq. (3) 
of Paper I, the peak surface brightness is 7^ ust = SjJ eak /Qbeam, 
where Qbeam is the solid angle of the telescope beam. As the dust 
temperature of the sources, we used the gas kinetic temperatures 
listed in Col. (4) of Table [1] and assumed that T^ust = 7\in- For 
the subcondensations, such as SMM 6b, 6c, etc., it was assumed 
that Tdust is the same as in the "main" core. This assumption 
may not be valid, however, because in most cases the subcon- 
densations lie outside the 40" beam of the NH3 measurements 
used to derive 7\i„ in the main cores. Moreover, the tempera- 
ture of an individual small condensation may be lower than the 
temperature of the parent core because of more effective shield- 
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Fig. 5. C 0(2 - 1) lines at 219 560.3568 MHz (JPL) arising from the image sideband. The x-axis range is different in the spectrum 
towards Ori B9 N, and the lines are overlaid with Gaussian fits for a better illustration. The blended OCS and CH3NH2 lines at 
~ 231.06 GHz marked on the spectrum towards Ori B9 N can be excluded on the basis of high upper-state energies and formation 
chemistry of the species (see text). The line is likely to be caused by an additional C ls O(2 - 1) velocity component. 



ing from the external radiation field. If the temperature is lower 
than assumed the mass and column density will be underesti- 
mated. As the dust opacity per unit dust mass at 350 /mi, #350 ^m, 
we used the value 1.0 m 2 kg -1 which was extrapolated from 
the Ossenkopf & Henning (1994, hereafter OH94) model de- 
scribing graphite-silicate dust grains that have coagulated and 
accreted thick ice mantles over a period of 10 5 yr at a gas den- 



sity of «h = "(H) + 2n(H2) - 2n(H2) = 10 5 crrT 3 . The same 
dust model was adopted in Paper I (/f870/«n — 0.17 m 2 kg -1 ), 
and is expected to be a reasonable model for cold, dense mole- 
cular cloud corefl We note that the sub-mm opacities recently 



7 For comparison, for dust grains covered by thin ice mantles at a 
density 10 5 cm 4 , £350^ - 0.78 m 2 kg -1 . The K 350/im value we have 
adopted is the same as in the OH94 model of grains with thin ice man- 
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Table 5. Spectral-line parameters. 



Source 


Transition 


V LSR 


Av 


7\lB 


J r MB dv° 


b i 
To ' 


7" ex c ' 








[km s~'] 


[km s~'[ 


[K] 


[Kkm s~'] 




[K] 


IRAS 05399-0121 


C"0(2- 


1) 


8.68 ± 0.02 


0.74 ± 0.06 


0.79 ±0.18 


1.27 ±0.14 [6.40, 9.64] 


0.10 ±0.05 


12.8 ± 1.5 




H 13 CO + (4 


-3) 


8.62 ± 0.08 


0.82 ±0.16 


0.09 ± 0.02 


0.08 ± 0.02 


0.11 ±0.05 


5.5 ±0.5 




DCO + (4 - 


•3) 


8.58 ± 0.03 


0.69 ± 0.06 


0.33 ± 0.04 


0.25 ± 0.03 


0.51 ± 0.10 


4.9 ±0.3 




N 2 H+(3 - 


2) 


8.73 ± 0.02 


0.71 ± 0.07 


0.88 ± 0.11 


1.00 ±0.11 [92.6%] 


2.80 ± 0.20 


5.0 ±0.2 




N 2 D+(3 - 


2) 


8.46 ± 0.04 


0.53 ± 0.02 


0.28 ± 0.03 


0.14 ±0.02 [91.6%] 


0.68 ±0.10 


5.0±0.2 rf 




o-D 2 CO(4 . 4 


-3 , 3 ) 


8.70 ± 0.22 


2.20 ± 0.46 


0.08 ± 0.03 


0.18 ±0.04 


0.01 ±0.001 


18.6 


SMM 1 


C 17 0(2 - 


1) 


9.20 ± 0.01 


0.93 ± 0.03 


0.96 ±0.17 


1.85 ±0.20 [7.00, 11.33] 


0.15 ±0.05 


11.3 ± 1.0 




H 13 CO + (4 


-3) 






< 0.13 










DCO + (4 - 


■3) 


9.31 ± 0.01 


0.47 ± 0.05 


0.55 ± 0.07 


0.42 ± 0.05 


1.89 ± 0.20 


4.5 ± 0.3 




N 2 H + (3 - 


2) 


9.31 ± 0.02 


0.66 ± 0.09 


0.65 ± 0.07 


0.70 ± 0.08 [92.6%] 


2.43 ± 0.10 


4.6 ± 0.3 




N 2 D + (3 - 


2) 


9.20 ± 0.03 


0.75 ± 0.07 


0.38 ± 0.04 


0.40 ± 0.04 [92.6%] 


1.71 ± 0.20 


4.6±0.3 rf 




o-D 2 CO(4 .4 


-3 0J ) 


9.21 ± 0.06 


0.67 ±0.16 


0.15 ± 0.02 


0.10 ±0.02 


0.02 ± 0.002 


18.6 


SMM 3 


C 17 0(2 - 


1) 


8.68 ± 0.06 


0.58 ± 0.09 


0.35 ± 0.05 


0.51 ± 0.09 [6.56, 9.52] 


0.05 ± 0.03 


11.0 ± 1.0 




DCO + (4- 


•3) 


8.54 ± 0.03 


0.42 ±0.16 


0.20 ±0.02 


0.09 ± 0.01 


0.24 ± 0.03 


5.1 ± 0.3 




N 2 H + (3 - 


2) 


8.57 ± 0.03 


0.85 ± 0.09 


0.62 ± 0.09 


0.67 ± 0.08 [92.6%] 


1.23 ± 0.10 


4.9 ± 0.3 




N 2 D+(3 - 


2) 


8.39 ± 0.04 


0.53 ±0.11 


0.21 ± 0.03 


0.17 ± 0.03 [91.6%] 


0.59 ± 0.10 


4.9 ± 0.3'' 


IRAS 05405-0117 


C 17 0(2 - 


1) 


9.20 ± 0.04 


0.54 ± 0.04 


0.30 ± 0.03 


0.17 ± 0.03 [78.7%] 


0.05 ± 0.03 


10.5 ± 0.5 


2nd v-comp. 


C"0(2 - 


1) 


1.30 ±0.33 


0.58 ± 1.08 


0.17 ± 0.02 


0.14 ±0.02 [78.7%] 






3rd v-comp. 


C"0(2 - 


1) 


2.98 ± 0.33 


0.54 ± 1.08 


0.20 ± 0.02 


0.13 ± 0.02 [78.7%] 








DCO + (4- 


•3) 






< 0.08 










N 2 H+(3 - 


2) 


9.30 ± 0.03 


0.69 ± 0.06 


0.66 ± 0.09 


0.56 ± 0.07 [90.9%] 


4.26 ±0.10 


4.5 ± 0.3 




N 2 D + (3 - 


2) 


8.93 ±0.15 


0.53 ±0.34 


0.06 ± 0.01 


0.04 ±0.03 [91.5%] 


0.16 ±0.02 


4.5 ±0.3 rf 


3rd v-comp. 


N 2 D + (3 - 


2) 


2.98 ± 0.20 


0.53 ± 0.56 


0.05 ± 0.01 


0.05 ± 0.02 [92.4%] 






SMM 4 (2nd v-comp.) 


C 17 0(2- 


1) 


1.67 ±0.05 


0.92 ±0.19 


0.32 ± 0.06 


0.61 ± 0.08 


0.06 ± 0.03 


9.7 ± 1.2 


2nd v-comp. 


H 13 CO + (4 


-3) 


1.53 ±0.08 


0.43 ±0.15 


0.15 ±0.02 


0.07 ± 0.02 


1.02 ±0.60 


4.0 ±0.2 


2nd v-comp. 


DCO + (4 - 


■3) 


1.65 ±0.03 


0.49 ± 0.09 


0.21 ± 0.04 


0.19 ±0.03 


0.76 ± 0.50 


4.0 ± 0.3 


2nd v-comp. 


N 2 H + (3 - 


2) 


1.62 ±0.78 


0.61 ±21.10 


0.32 ± 0.07 


0.27 ± 0.07 [92.6%] 


1.59 ±0.50 


4.0 ±0.3 


2nd v-comp. 


N 2 D + (3 - 


2) 


1.49 ±0.04 


0.54 ±0.10 


0.24 ± 0.03 


0.17 ±0.03 [90.9%] 


1 .50 ± 0.20 


4.0 ± 0.3 rf 


SMM 5 


C 17 0(2 - 


1) 


9.37 ± 0.03 


0.54 ±0.01 


0.22 ± 0.02 


0.14 ±0.02 [78.7%] 


0.04 ± 0.03 


10.1 ±0.5 




DCO + (4 - 


■3) 






< 0.07 










N 2 H+(3 - 


2) 


9.43 ± 0.07 


0.44 ±0.17 


0.16 ±0.03 


0.14 ±0.02 [92.6%] 


0.67 ± 0.04 


3.8 ±0.1 




N 2 D+(3 - 


2) 






< 0.11 








SMM 6 


C 17 0(2 - 


1) 


9.51 ±0.04 


0.54 ± 0.04 


0.28 ± 0.03 


0.15 ±0.03 [69.4%] 


0.05 ± 0.03 


9.7 ± 0.5 




DCO + (4 - 


•3) 


9.47 ± 0.02 


0.42 ± 0.04 


0.37 ± 0.05 


0.19 ± 0.02 


3.47 ± 0.20 


4.0 ±0.3 




N 2 H+(3 - 


2) 


9.52 ±0.03 


0.65 ± 0.09 


0.51 ±0.06 


0.44 ± 0.06 [92.6%] 


5.75 ±0.35 


4.2 ± 0.2 




N 2 D + (3 - 


2) 


9.25 ± 0.03 


0.56 ±0.10 


0.33 ± 0.03 


0.26 ± 0.03 [90.9%] 


2.92 ± 0.30 


4.2±0.2 rf 


OriB9N 


C 17 0(2 - 


1) 


9.29 ±0.33 


0.54 ± 1.08 


0.30 ± 0.05 


0.20 ± 0.03 [96.0%] 


0.04 ± 0.03 


11.6 ± 1.0 


2nd v-comp. 


C 17 0(2 - 


1) 


1.85 ±0.33 


1.39 ± 1.08 


0.29 ± 0.05 


0.70 ± 0.08 


0.04 ± 0.03 


11.7 ±2.0 




DCO + (4 - 


■3) 


9.42 ± 0.19 


0.47 ± 1.38 


0.04 ± 0.01 


0.02 ± 0.02 


0.05 ± 0.02 


4.9 ± 0.3 


2nd v-comp. 


DCO + (4 - 


■3) 


2.11 ±0.11 


0.78 ±0.31 


0.09 ± 0.01 


0.09 ± 0.02 


0.14 ±0.08 


4.6 ± 0.5 




N 2 H+(3 - 


2) 


8.80 ±0.16 


0.44 ± 0.24 


0.05 ± 0.02 


0.05 ± 0.02 [92.6%] 


0.12 ±0.03 


4.1 ±0.2 


2nd v-comp. 


N 2 H+(3 - 


2) 


1.81 ±0.11 


0.85 ±0.11 


0.09 ± 0.02 


0.13 ±0.03 [92.6%] 


0.23 ± 0.07 


4.1 ±0.2 




N 2 D + (3 - 


2) 






< 0.07 








SMM 7 


C 17 0(2 - 


1) 


3.75 ± 0.02 


1.04 ±0.05 


0.66 ± 0.09 


0.82 ±0.09 [1.38, 4.95] 


0.15 ±0.05 


8.9 ± 1.0 


"9 km s _1 "-comp. 


C 17 0(2 - 


1) 


9.04 ± 0.07 


0.55 ± 1.74 


0.10 ±0.06 


0.07 ± 0.03 [69.4%] 








DCO + (4 - 


•3) 


3.70 ± 0.04 


0.42 ±0.10 


0.22 ± 0.02 


0.10 ±0.02 


0.94 ± 0.20 


3.9 ± 0.3 




N 2 H+(3 - 


2) 


4.01 ±0.14 


0.67 ± 0.27 


0.14 ±0.02 


0.17 ±0.04 [92.6%] 


0.49 ± 0.07 


3.9 ±0.2 




N 2 D+(3 - 


2) 






< 0.07 









Notes. <n) Integrated intensity is derived from a Gaussian fit or, in the case of some C 17 lines, by integrating over the velocity range indicated in 
brackets. The percentage in brackets indicates the contribution of hf component's intensity lying within the Gaussian fit.'*' For C n O, H 13 CO + , 
DCO + , N 2 H + , and N 2 D + r is the optical thickness in the centre of a hypothetical unsplit line (see text). (c) For C 17 0, H I3 CO + , DCO + , and 
N2H + the values of To and T cx were estimated using RADEX with the kinetic temperature (7\j„) and the H2 density ((;i(H 2 )». The value of To for 
N 2 D + and o-D 2 CO was estimated using CLASS/Weeds; see Sect. 4.4 for details. {d) T cx is assumed to be the same as for N 2 H + (3 - 2). 



calculated by Ormel et al. (2011) are comparable to the OH94 
values (after 10 5 yr of coagulation). However, the dust opaci- 
ties are likely to be uncertain by a factor of > 2 (e.g., OH94; 
IMotte & Andre 20011 IQrmel etal. 201 U . For the average dust- 
to-gas mass ratio, = (M,i ust /M gSL! ), we adopted the canonical 
value 1/100. Finally, we assumed a He/H abundance ratio of 0.1, 
which leads to the mean molecular weight per H2 molecule of 

m 2 = 2.8. 

The integrated flux densities used to calculate the masses 
refer to core areas, A. Thus, in order to properly calculate the 
volume-average H2 number density, (n(H2)), we use the effec- 
tive radius, ^ e ff = y/A/n, in the formula 

<n(H 2 )> = J^-, (i) 

m 2 m n 



ties at a density 10 6 cm 3 . At 10 6 cm 3 with thick ice mantles k 350 
rises to about 1.1m 2 kg -1 . 



where (p) = M/ ^47r/3 x /r ff j is the mass density, and mh is the 
mass of a hydrogen atom. 

The results of the above calculations are presented in Table|6] 
The uncertainties in the derived parameters were propagated 
from the uncertainties in 7\m- We note that the lcr rms noise on 
our SABOCA map, ~ 0.06 Jy beam -1 , corresponds to a 3<x mass 
detection limit of ~ 0.1 Mq assuming 7d U st = 10 K. In terms of 
column density, the 3<x detection limit is about 2.8 x 10 21 cm -2 . 



4.2. Dust properties determined from 350 and 870 /jm data 

We can use our observations at two submm wavelengths, 350 
and 870 /mi, to estimate the dust colour temperature, r^st, and 
dust emissivity spectral index, /3 [/q = Ko(Ao/A~f]. Note that ac- 
cording to the Wien displacement law, the wavelength of the 
peak of the blackbody radiation curve is about 290 //m at 10 K. 
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IRAS 05399-0121 
D 2 CO(4 04 -3 OJ ) 
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1 J 
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v LSR ftm s" 1 ] 

SMM 1 

D Z CO(4 0|4 -3 OJ ) . 



V 



v«n [km • J 

Fig. 4. orf/!o-D 2 CO(4 , 4 - 3 ,3) spectra towards IRAS05399 and 
SMM 1. 



Table 6. The effective radius, mass, H2 column density, and 
volume-averaged H2 density derived from the 350 fim emission. 


Source 


Keff 


M 


N(R 2 ) 


<«(H 2 )) 




[pc] 


[Mq] 


[10 22 crrT 2 ] 


[10 5 crrT 3 ] 


SMM 3 


0.03 


2.1 ±0.8 


10.4 ± 2.7 


4.0 ± 1.5 


SMM 3b 






0.7 ± 0.2" 




SMM 3c 






0.7 ± 0.2" 




IRAS 05405-0117 


0.02 


0.2 ±0.1 


1.3 ±0.3 


2.0 ± 1.0 


SMM 4 


0.01 


0.04 ± 0.02 


0.4 ±0.1 


0.9 ± 0.5 


SMM 4b 


0.02 


0.1 ±0.04° 


0.7 ±0.1" 


1.0 ±0.4 


SMM 5 






0.6 ±0.1 




SMM 6 


0.02 


0.2 ±0.1 


1.0 ±0.1 


1.6 ±0.8 


SMM 6b 


0.02 


0.1 ±0.05" 


0.9 ±0.1" 


1.3 ±0.6 


SMM 6c 


0.01* 


0.1 ±0.04"-* 


0.7 ±0.1" 


1.9 ±0.7* 


SMM 6d 






0.7 ±0.1" 




Ori B9 N 






0.3 ±0.1 




SMM 7 


0.02 


0.6 ±0.3 


1.9 ± 1.0 


2.6 ± 1.3 


SMM 7b 






1.3 ±0.7" 





Notes. (fl) Calculated by assuming T imt is the same as for the "main" 
core.'*' These values include contributions from both SMM 6b and 6c. 



Therefore, the Rayleigh-Jeans (R-J) approximation, hv <sz k^T, 
is not valid for our 350 and 870 //m data. 

In the following analysis, it is assumed that T^st and ft are 
constant across the source. We first smoothed the SABOCA 
map to the resolution of our LABOCA data, and then calculated 
the flux densities at both wavelengths in a fixed 40" diameter 
aperture [Col. (7) of Table [3) . The resulting flux density ratios, 

^350 /^870 ' are §i ven m Col. (2) of Table [7] A dust colour tem- 
perature can be determined by fixing the value of ft [see, e.g., 
Eq. (3) of Shetty et al. (2009)]. The value of ft in the OH94 thick- 
ice dust model we adopted earlier is about 1 .9 over the wave- 
length range A e [250, 1300 /an] (see also |Shirley etal. 2005| ). 
By adopting the value ft — 1.9, we derive the T& ast values in the 
range ~ 7.9 - 10.8 K [see Col. (3) of TableQ. 

To calculate ft from the ratio of two flux densities at different 
wavelengths, an estimate for the dust temperature is needed. We 
assumed that Tjust = 7kim and applied Eq. (5) of Shetty et al. 
(2009). The resulting values, ft » 0.5 - 1.8, are listed in Col. (4) 
of Table [7] We note that T^\ n measurements were obtained with 
a 40" resolution, whereas the above flux density ratio was de- 
termined at about 20" resolution. However, flux densities were 
measured in a 40" aperture, which matches the resolution of our 
NH 3 data. 



Table 7. Dust colour temperature and emissivity index derived 
from 350-to-870 //m flux density ratio. 



Source 


J ISfl 


,r40" 
' J 87(1 


Tdust" [K] 




SMM 3 


2.4 


± 1.1 


10.8^6 


1.8 ±0.6 


IRAS 05405-0117 


1.5 


±0.6 


8.7 


1.1 ±0.6 


SMM 4 


1.4 


+ 0.8 


8.4 +i5 


0.5 ±0.8 


SMM 5 


1.7 


±0.9 


9.0+11 


1.2 ±0.8 


SMM 6 


1.2 


±0.6 


7 - 9+ is 


0.8 ±0.9 


Ori B9 N 


1.7 


±0.9 


9.0+" 
8 9+^ 


0.8 ±0.8 


SMM 7 


1.6 


±0.9 


1.7 ±0.8 



Notes. <a) Derived by assuming ft = 1.9 (see text).' 1 Calculated by 
assuming T iust = T Un as derived from NH 3 observations. 

4.3. SEDs 

With the aid of our new 350 /mi data, we were able to refine some 
of the protostellar SEDs presented in Paper I. As in Paper I, the 
SEDs constructed from the 24, 70, 350, and 870-//m flux den- 
sities were fitted by the sum of two modified blackbody curves 
of cold and warm temperatures (the original version of the fit- 
ting routine was written by J. Steinacker). Again, the dust model 
used to fit the SEDs is the OH94 model of thick ice mantles as 
adopted earlier in Sect. 4.1, incorporating a dust-to-gas mass ra- 
tio of 1/100. The SEDs of SMM 3, IRAS05405, and SMM 4 
are shown in Fig. [6] Note that in Paper I we utilised the IRAS 
flux densities to build the SED of IRAS05405. In the present 
paper we have ignored the IRAS data from the fit because the 
~arcminute resolution of the IRAS observations also includes 
emission from the nearby core SMM 4, confusing the emission 
from IRAS05405. This probably explains why the 100-/mi flux 
density of IRAS05405 appears so high (~ 19.7 Jy). 

The SED fitting results are shown in Table [8] Column (2) of 
Table[8]gives the mass of the cold component, which represents 
the mass of the cold envelope, M co id = M env . In Cols. (3) - (5), 
we list the luminosities of the cold and warm component, and the 
bolometric luminosity, Lboi = ^coid + ^wai-m- The temperature of 
the cold component, r co id, is given in Col. (6). Columns (7) and 
(8) give the Lcoid/^boi and L S ubmm/£boi ratios, where the submm 
luminosity, L su b mm , is defined to be the luminosity longward of 
350 jt/m. In the last column of Table [8] we give the normalised 
envelope mass, M env /L^ (Bontemps et al. 1996). The latter pa- 
rameter, which is related to the outflow activity, decreases with 
time and can be used to further constrain the core evolutionary 
stage. We do not report the temperature and mass of the warm 
component, because i) the OH94 dust model of grains covered 
by thick ice mantles we have adopted may not be appropriate for 
the warm dust component; and ii) dust emission may not be opti- 
cally thin at 24 and 70 /mi, so the masses and temperatures of the 
warm component are not well constrained. In contrast, the lumi- 
nosity of the warm component is not affected by these opacity 
effects. 

4.4. Molecular column densities and fractional abundances 

The beam-averaged column densities of C n O, H 13 CO + , DCO + , 
and N2H + were derived using a one-dimensional spherically 
symmetric non-LTE radiative transfer code called RADEX (see 
Sect. 3.3). RADEX uses the method of mean escape probability 
for an isothermal and homogeneous medium. The molecular 
data files (collisional rates) used in the RADEX excitation analy- 
sis were taken from the LAMDA database ( Scho ier et al. 2005l l. 
The C n O, H 13 CO + , DCO + , and N 2 H + transitions are treated as 
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1 10 100 1000 10000 1 10 100 1000 10000 1 10 100 1000 10000 

A [/Am] A [M 111 ] A [/Am] 



Fig. 6. Spectral energy distributions of three protostellar cores in Orion B9 built from Spitzer 24 and 70 pm, S ABOCA 350 pm, and 
LABOCA 810-pm flux densities. Open circles in the middle panel represent IRAS data points at 12, 25, 60, and 100 /mi (not used 
in the fit). The solid lines correspond to the two-temperature model fits to the data. Error bars (lcr) are indicated for all data points, 
but are mostly smaller than the symbol size. Note the appearance of a 10-p.m silicate absorption feature in the SED of IRAS05405, 
and the absorption "knee" at ~ 30 pm (most notably towards IRAS05405) where a considerable change of k v occurs at the dust 
model used (OH94 and references therein). 



Table 8. Fitting results of the SEDs. 



Source 




^cold 


^ warm 


r h 

M>ol 




icold/^bol 


^submm / ^bol 


M 1 1 °- 6 




[M ] 


[Lq] 


[Lq] 


[L©] 


[K] 






[M Q /L£f] 


SMM3 


8.2 ±2.6 


0.3 ±0.1 


0.9 ±0.1 


1.2 ±0.1 


8.0 


0.3 ±0.1 


0.1 


7.4 ±3.2 


IRAS 05405-0117 


1.6 + 0.5 


0.06 ± 0.02 


2.3 ±0.2 


2.3 ±0.2 


8.0 


0.03 ± 0.01 


0.01 


1.0 ±0.3 


SMM4 


2.0 ±0.8 


0.07 ± 0.03 


0.5 ± 0.02 


0.6 ± 0.04 


8.0 


0.1 ±0.05 


0.05 


2.7 ± 1.1 



Notes. w M cold = M env . (i,) L bo i = L warm + L cM . 



a hypothetical unsplit transition. The input parameters in the off- 
line mode of RADEX are the gas kinetic temperature, H2 num- 
ber density, and the width (FWHM) and intensity of the spectral 
line. We used the values of 7\j n and (n(H2)) listed in Table [1] 
However, we multiplied the densities by 1.2 (He/H2 = 0.2) to 
take the collisions with He into account (see Sect. 4.1 of the 
RADEX manuafl). As the input line intensity we used the main- 
beam brightness temperature, Tmb- When the source is resolved, 
7mb is equal to the R-J equivalent radiation temperature, Tr. 
The simulations aim to reproduce the observed line intensity 
and yield the values of To and T ex , and the total column den- 
sity of the molecule (N tot ). We varied 7\j n and (n(H2)) according 
to their errors to estimate the uncertainties associated with to, 
r ex , and N tot . The uncertainties in Tmb and Av were not taken 
into account, but test calculations showed that they lead to er- 
rors in column density that are comparable to those derived from 
the errors in temperature and density. We also note that the use 
of higher H2 densities for the cores derived from the SABOCA 
map [Col. (4), Table |6l would lead to lower column densities 
of the molecules because then the excitation would be closer to 
therm alisation. 

For N2D + and 0-D2CO there are no molecular data files 
available in the LAMDA database. The line optical thicknesses 
and total beam-averaged column densities of these molecules 
were determined through LTE modelling with CLASS/Weeds. 
The input parameters for a Weeds model are N to t, Tex, source size 
(0 S ), linewidth (FWHM), and offset from the reference-channel 
velocity. The linewidth is directly determined from the observed 
line profile, so there are basically three free parameters left (Mot, 
T ex , S ). Some of the model parameters may be degenerate, 
and cannot be determined independently (Sc hilke et al. 2006; 
IMaret et al. 201 II ). The source size is degenerate with excita- 



tion temperature in the case of completely optically thick lines 
(r » 1), and with column density if the lines are completely 
optically thin (r <K 1). We assumed that the source fills the tele- 
scope beam, i.e., that the beam filling factor is unit)Q, and thus 
the line brightness temperature is Tb - 7mb [see Eqs. (1) and 
(2) in Maret et al. (201 1)]. For N2D + , we used as T ex the values 
obtained for N2H + from RADEX simulations. For the asymmetric 
top rotor 0-D2CO we adopted the value T ex = 2/3x__ u /£b, which 
gives a lower limit to the column density ( Hate hell et al. 1 998). 
The input N tot was then varied until a reasonable fit to the line 
was obtained (see Fig - [7j . The associated error estimate is based 
on the 10% calibration uncertainty. In the derivation of column 
densities, the line-strength contribution of hf components within 
the detected lines was taken into account, i.e., the column densi- 
ties were corrected for the fraction of the total line strength given 
in brackets in Col. (6) of Table[5] 

We also determined the HCO + column density from the 
column density of H 13 CO + . For this calculation, it was as- 
sumed that the carbon-isotope ratio is [ 12 C]/[ 13 C] = 60 
dWilson & Rood 19941 [Savage et al. 2002| ). This value has been 
used in several previous studies, e.g., by Bergin et al. (1999) in 
their study of dense cores in Orion (including IRAS05399C3- 



9 The 27" beam of our N 2 D + observations is comparable to the ex- 
tent of the strongest dust emission region within the cores. On the other 
hand, N 2 D + emission has been found to trace dust emission very well in 
low-mass dense cores (e.g., Cra psi et al. 2005| l. Therefore, the assump- 
tion of unity beam filling factor is reasonable. For o-D 2 CO, the beam 
filling factor may be < 1. For example, Bergman et al. (201 1), using the 
same resolution as we (27"), derived the filling factor of = 0.5 for the 
o-D2CO(4o.4 - 3o.3) emission region in p Oph A. Moreover, the filling 
factor is lower if the cores contain unresolved small-scale structure. 

10 Spectral lines of 12 C-isotopologue of HCO + are likely to be opti- 
cally thick. Therefore, the HCO + deuteration can be better investigated 

http://www.sron.rug.nl/~vdtak/radex/radex_manual.pdf through the DCO + /H 13 CO + column density ratio. However, a caveat 
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We calculated the fractional abundances of the molecules 
by dividing the molecular column density by the H2 column 
density: x(mol) = N(mol)/N(H2). For this purpose, the va- 
lues of ./V(H2) were derived from the LABOCA dust continuum 
map smoothed to the corresponding resolution of the line ob- 
servations. The resolution of the H 13 CO + observations (18") is 
slightly better, but comparable, to that of the original LABOCA 
data, and thus no smoothing was done in this case. The de- 
rived column densities and abundances are listed in Table [9] We 
stress that the reported uncertainties are formal and optimistic, 
and probably underestimate the true uncertainties. A stock chart 
showing the fractional abundances (excluding the additional ve- 
locity components) is presented in Fig. [8] The abundance errors 
were derived by propagating the errors in N(mol) and /V(H2). 



SMM 1 
N 2 D + (3-?) 


; 

\A ,aAA r 





Fig. 7. Example of the Weeds LTE modelling outlined in 
Sect. 4.4. The synthetic model spectrum is overlaid as a red line. 
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Fig. 8. A bar representation of fractional abundances on a log- 
arithmic scale. The H 13 CO + and HCO + abundances derived to- 
wards IRAS05399 are shown by a slightly stretch stock for a 
better illustration. 



4.5. CO depletion and deuterium fractionation 

To estimate the amount of CO depletion in the cores, we calcu- 
lated the CO depletion factor, /b- If ■x(CO) can is the "canonical" 



should be noted here. The HCO + molecules are produced directly from 
CO (Sect. 4.6). On the other hand, at low temperature, CO is suscepti- 
ble to the exothermic isotopic charge exchange reaction 13 C + + 12 CO — » 
12 C + + 13 CO + AE, where AE/k B = 35 K (Wats on et al.T976l . This is 
expected to cause considerable 13 C-fractionation in cold and dense gas, 
which complicates the deuteration analysis. 



(undepleted) abundance, and x(CO) b s is the observed CO abun- 
dance, fo is given by 



jc(CO) c 

4CO)„ 



(2) 



We adopted the standard value jc(CO) can = 
(Frerking et al. 1982). To calculate the "canonical 



9.5 x 1(T 5 
C 17 abun- 
dance we assumed the oxygen-isotopic ratio of [ 16 0]/[ 18 0] = 
500 (e.g., I Wilson & Rood 19941 IWilliams etal. 19981 1. When 
this is combined with the [ 18 0]/[ 17 0] ratio, for which we 



use the standard value 3.52 ( Frerking et al. 1982 1, the value of 
x(C n O) can can be calculated as 



*(C 17 0) can = 



i(CO) c 



[ 18 0]/[ 17 0] x [ 16 0]/[ 18 0] 



-*(CO) ca 
1760 



(3) 



The depletion factor fu is then calculated from /d = 
x(C 17 0) can /x(C 17 0) obs . The results are listed in Col. (2) of 
Table [TO] where the +-errors quoted were calculated by pro- 
pagating the uncertainty in x(C I7 0) b s . 

The degree of deuterium fractionation in HCO + and 
N2FF was calculated by dividing the column density of the 
deuterated isotopologue by its normal hydrogen-bearing form 
as fl D (HCO + ) ee jV(DCO + )//V(HCO + ) and fi D (N 2 H + ) = 
A^(N 2 D + )/A^(N 2 H + ). The error in R D was derived from the er- 
rors in the corresponding column densities [see Cols. (3) and (4) 
of Table [10|. We note that the N 2 H + and N 2 D + column densi- 
ties were calculated using the non-LTE and LTE models, respec- 
tively. Although not taken into account here, this may introduce 
an additional error of a factor of a few in Ro, which should be 
bear in mind. Note, however, that in all cases the normal and 
deuterated forms of the molecule show line emission at similar 
radial velocities and with similar linewidths. Therefore, the two 
transitions are probably tracing the same gas which makes the 
derived deuteration levels reasonable. 



4.6. Fractional ionisation 

Fractional abundances of all the observed ionic species and 
their different isotopologues could be determined only towards 
IRAS05399 and SMM 4-LVC. Therefore, the fractional ionisa- 
tion and cosmic -ray ionisation rate of H2 were estimated only 
for these two sources. A rough estimate of the lower limit to 
the ionisation degree can be obtained by simply summing up the 
abundances of the ionic species (e.g., Casel li et al. 2002al Paper 
I): 

x(e) > ^(HCO + ) + 4H 13 CO + ) + jc(DCO + ) + x(N 2 H + ) + a;(N2D + ). 

(4) 

This is based on the gas quasi-neutrality: the electron abun- 
dance equals the difference between the total abundances of the 
cations and anions. The resulting values are x(e) > 1 .5 x 10~ 8 for 
IRAS05399, and x(e) > 6.1 x 10" 8 for SMM 4-LVC. In the outer 
envelope where CO is not heavily depleted, HCO + is expected 
to be the main molecular ion. On the other hand, we do not have 
observational constraints on the abundances of H + , (and its 
deuterated isotopologues), H30 + , and metal ions (such as C + ), 
all of which could play an important role in the ionisation level. 

When the fractional ionisation in the source is determined, 
the abundance ratio Ru = [HCO + ]/[CO] can be used to infer 
the cosmic-ray ionisation rate of H2, £h 2 • By deriving a steady- 
state equation for the H3 abundance, and applying it in the cor- 
responding equation for HCO + , it can be shown that 
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Table 9. Molecular column densities and fractional abundances with respect to H2. 



Source 


2V(C lv O) 


N(N 2 H + ) 




N(H li CO + ) 


N(HCO + ) 


N(DCO + ) 


N(o - D 2 CO)° 




[10 14 cm" 2 ] 


[10 13 cm" 2 ] 


[10 12 cm- 2 ] 


[10 12 cm- 2 ] 


[10 14 cm- 2 ] 


[10 12 cm- 2 ] 


[10 12 cm" 2 ] 


IRAS 05399-0121 


3.4 ±0.5 


1.5 ±0.3 


3.1 ±0.3 


6.4 ± 1.0 


3.8 ±0.6 


7.6 ±0.9 


2.0 ±0.2 


SMM 1 


5.6 ± 0.5 


1.2 ±0.3 


11.9 + 1.2 






13.1 ± 1.0 


1.0 ±0.1 


SMM3 


1.3 ±0.3 


0.8 ±0.2 


2.7 ±0.3 






2.3 ±0.3 




IRAS 05405-0117 


1.3 ±0.4 


2.3 ±0.3 


0.8 ±0.1 










SMM 4 
















2nd v-comp.* 


2.0 ±0.5 


1.0 ± 0.3 


9.5 ± 1.0 


18.2 ± 7.8 


10.9 ± 4.7 


10.4 ± 5.0 




SMM 5 


1.0 ±0.4 


0.4 ±0.1 












SMM 6 


1.4 ±0.4 


3.0 ±0.3 


17.6 ± 1.8 






25.9 ±2.0 




Ori B9 N 


0.9 ±0.3 


0.1 c 








1.8 ±0.3 




2nd v-comp. 6 


2.3 ±0.5 


0.4 ±0.1 








6.4 ± 2.0 




SMM 7 


5.2 ±0.8 


0.5 ±0.1 








10.0 ±2.0 






x(C"0) 


x(NoFT) 


x(N 2 D + ) 


x(H 1J CO + ) 


x(HCO + ) 


x(DCO + ) 


x(o - Do CO) 




[10~ 8 ] 


[10- 10 ] 


[10- 10 ] 


[10- 10 ] 


[10~ 8 ] 


[10- 10 ] 


[io- H ] 


IRAS 05399-0121 


1.7 ±0.4 


6.5 ± 1.9 


1.5 ±0.3 


2.3 ± 0.6 


1.4 ±0.4 


3.6 ±0.9 


9.8 ±2.2 


SMM 1 


2.9 ± 0.5 


5.4 ± 1.6 


6.1 ± 1.1 






5.8 ±0.9 


5.1 ±0.9 


SMM 3 


0.5 ±0.1 


3.5 ± 1.0 


1.1 ±0.2 






1.0 ±0.2 




IRAS 05405-0117 


1.5 ±0.5 


23.9 ±4.0 


0.9 ±0.2 










SMM 4 
















2nd v-comp.* 


1.5 ±0.5 


6.2 ±2.5 


6.8 ± 2.0 


9.6 ±4.8 


5.8 ±2.9 


6.3 ±3.5 




SMM 5 


1.5 ±0.6 


5.5 ± 1.5 












SMM 6 


1.3 ±0.4 


22.3 ± 2.7 


15.4 ± 1.9 






18.7 ± 1.9 




Ori B9 N 


2.6 ± 1.0 


3.3 ±0.6 








6.3 ± 1.5 




2nd v-comp.* 


6.7 ± 2.6 


13.4 ±5.4 








22.9 ± 10.2 




SMM 7 


3.3 ± 1.0 


2.8 ±0.9 








5.5 ± 1.8 





Notes. <a) The N 2 D + and o-D 2 CO column densities were derived from Weeds models. For o-D 2 CO we assumed that T a = 2/3 x E n lk^. ib) For 
SMM 4-LVC (~ 1.6 km s" 1 ) and Ori B9 N-LVC (~ 1.9 km s" 1 ), T kin = 10.4 ± 1.4 K and 13.6 ± 2.5 K, respectively (Paper II). <c) The associated 
error derived by varying 7\i n and («(H 2 )) in the RADEX calculation is negligible. 



[fa 2 /«(H 2 )]^ 

[pn+x(e) + k H +x{CO) + £ gl , H +x(g)] [y8 H co+x(e) + £ sr , H co+x(g)] 

(5) 

where ku+ is the rate coefficient for the reaction + 

CO — > HCO + + H2, j8h+ and Pnco+ are the dissociative re- 
combination rate coefficients of H3 and HCO + , and £ gr H+ and 
fegr,HCO + are the rate coefficients for the recombination of and 
HCO + onto dust grains. The values of ku+ and /3hco+ were taken 
from the UMIST databaseO dWoodall et al. 20071 1. whereas j8 H +, 
^gr,H + , and ^ gr> HCO+ were interpolated from Pagani et al. (2009a; 
Tables A.l and B.l therein). For the grain abundance, x(g), we 
used the value 2.64 x lO -12 , which is based on the grain radius 
a = 0.1 //m, density p gra i n = 3 g cirr 3 , and the dust-to-gas mass 
ratio R d = 1/100 [see, e.g., Eq. (15) in Pagani et al. (2009a)]. 
To calculate we adopted as x(e) the summed abundance of 
ionic species. As MHt) of SMM 4-LVC, we used the density of 
SMM 4. The obtained values are ~ 2.6 x 10" 17 s -1 towards 
IRAS05399, and ~ 4.8 x lO" 16 s" 1 towards SMM 4-LVC. These 
values should be taken as lower limits in the sense that we have 
used the lower limits to x(e). 

5. Discussion 

5.1. Dust properties 

By fixing the value of ft to 1 .9, we derived the dust temperatures 
in the range ~ 7.9- 10.8 K. These are about 0.5-5.5 K lower than 
the gas temperatures in the same objects as derived from NH3. 

11 http://www.udfa.net/ 



Table 10. The degree of CO depletion and deuterium fractiona- 
tion. 



Source 


fi> 




^ D (HCO + ) 


fl D (N 2 H + ) 


IRAS 05399-0121 


3.2 ± 


0.7 


0.020 ± 0.004 


0.207 ± 0.046 


SMM 1 


1.9 ± 


0.3 




0.992 ± 0.267 


SMM 3 


10.8 i 


.2.2 




0.338 ± 0.092 


IRAS 05405-0117 


3.6 ± 


1.2 




0.035 ± 0.006 


SMM 4 










2nd v-comp. 


3.6 ± 


1.2 


0.010 ± 0.006 


0.950 ± 0.302 


SMM 5 


3.6 ± 


1.4 






SMM 6 


4.2 ± 


1.3 




0.587 ± 0.084 


Ori B9 N 


2.1 ± 


0.8 






2nd v-comp. 


0.8 ± 


0.3 






SMM 7 


1.6 ± 


0.5 







Due to the large uncertainties associated with Td asl , the near 
equality 7\j n - 7d us t seems possible, as expected at high den- 
sities where collisional coupling between the gas and dust be- 
comes efficient. Theoretical models have shown that in the dense 
interiors of starless cores [niWi) > 3 x 10 4 cm- 3 ], the gas and 
dust temperatures are similar, although the gas can be slightly 
warmer due to cosmic -ray heating (e.g. JGaili et al. 20021 . 

Some of the dust temperatures we derived are very low. 
In particular, for SMM 6 we obtained r dust ^ 7.9+ 3 g K. 
Theoretical models (e.g., Evans et al. 2001) and previous ob- 
servational studies (e.g. JSchnee et al. 2007al |Crapsi et aL~2 007; 
Harju et al. 2008) > have indicated that very low gas and dust tem- 
peratures (~ 6 - 7 K) can be reached in the dense interiors of 
dense cores. 

The dust emissivity spectral indices we derived, (3 ~ 
0.5 - 1.8, are physically reasonable, and suggest that the as- 
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sumption T^t = 7\in is valid (cf. ISchnee & Go odman 2005 ). 
Furthermore, the ft values derived towards SMM 3, 5, and 7 are 
close to 1 .9. We note that the dense gas and dust associated with 
the additional velocity components along the line of sight also 
affect the observed dust continuum properties (e.g., SMM 4). 
Therefore, the derived ft values for these sources should be taken 
with caution. 

We note, however, that a decrease of ft from the 'fiducial' 
value 2 to a shallower emissivity spectral index could be re- 
lated to dust grain coagulation in the inner parts of dense cores 
(e.g., Miyake & Nakagawa 1993; OH94). Observational studies 
have found evidence that ft decreases as a result of grain growth 
at high densities (e.g., IGoldsmith et al. 1997IIVisser et al. 1998b . 
More recently, Kwon et al. (2009) found that ft < 1 for their 
sample of three Class sources, resembling the values we found 
towards the Class candidates IRAS05405 and SMM 4. These 
results suggest that dust grains in the envelopes of Class pro- 
tostars can grow in size leading to a shallower spectral index 
of dust emissivity. Based on Herschel observations of cold in- 
terstellar clouds detected with the Planck satellite, Juvela et 
al. (2011) found that ft decreases down to ~ 1 near internal 
heating sources. Radiative transfer modelling, however, sug- 
gests that such a decrement is due to line-of-sight tempera- 
ture variations rather than changes in the grain properties. On 
the other hand, some studies of low-mass dense cores sug- 
gest that ft could be larger than 2 (e.g., Shirley et al. 2005 



ISchnee et al. 2010"al Shirley et al. 201 1 1. Also, recent studies of 
Planck-detected cold Galactic clumps indicate that ft > 2 
dPlanck Co llaborati on et al. 201 11 . It is possible that ft is anti- 
correlated with 7d ust , so that in cold, dense cores the emissivity 
spectral index is steeper than 2. 



5.2. Refined SEDs 

Inclusion of the new 350 fim data to the source SEDs con- 
firmed our earlier protostellar classifications. In Paper I, we 
classified the sources SMM 3, SMM 4, and IRAS05405 as 
Class objects. For such objects, the bolometric temperature is 
Thai < 70 K, L S ubmm/£boi ratio is > 0.005, and they are charac- 
terised by the ratio M env /L° o 6 j > 0.4 M©/L^ 6 dAndreetal. 19931 



Bontemps et al. 1996 I Andre et al. 20001 1. All these conditions 
are fulfilled after the refined SED analysis. 

The SED models presented in Paper I suggested 350-^m flux 
densities of 21.7 Jy for SMM 3, 13.5 Jy for IRAS05405, and 
1 1 . 1 Jy for SMM 4 - much higher than determined in the present 
work. Compared to the SED results in Paper I, the envelope mass 
of SMM 3 and 4 is 1.0 M© higher and 1 .8 M© lower, respec- 
tively; the mass of IRAS05405 remains the same. The refined 
luminosity is lower by a factor of ~ 3 in all three cases. The 
temperature of the envelope, r co id, resulting from the new SEDs 
is only 8 K for all sources. This is much lower than the values 
1 1 .6-16. 1 K derived in Paper I. The L C oid/^boi ratios derived here 
are also much smaller than estimated previously. Unlike deduced 
in Paper I, the cold component does not appear to dominate the 
total source luminosity. The L su \, mm /L\, i and M em /LP^ ratios de- 
rived here are mostly comparable to those obtained in Paper I. 

The angular resolution of the Spitzer images used is about 
6" at 24 jt/m, and 18" at 70 //m. Photometry was done by point- 
source fitting, and the aperture size for photometry was 13" at 
24 yum, and 35" at 70 fj.m (Paper I). The submm data points used 
in the SEDs were obtained at about 20" resolution within a 40" 
aperture. Therefore, the 70, 350, and 870-/im flux densities refer 
to a similar spatial scale (~ 0.09 pc). Direct comparison of these 



data with the 24 fim data is reasonable because we are dealing 
with 24-yizm point sources. Our model SEDs therefore produce an 
approximation to the core parameters on ~ 0.09 pc spatial scale. 
Note that the sources SMM 4 and 4b are treated as one source 
and thus the corresponding SED should be taken with caution 
(Sect. 5.7.2). 

5.3. Molecular column densities and abundances 

The present N2I-F column densities are about 0.6-5 times the 
corresponding LTE column densities presented in Paper II; 
within the errors, the two values are comparable in the case of 
IRAS05399, SMM 1, 3, and 7. Bergin et al. (1999) found the co- 
lumn densities JV(N 2 H + ) ~ 4.0 + 0.3 x 10 12 cnr 2 , JV(H 13 CO + ) ~ 
6.9 ± 1.3 x 10 11 cm" 2 , and JV(DCO + ) ~ 1.4 ± 0.2 x 10 12 cnr 2 
towards IRAS05399 using a statistical equilibrium model. Our 
values, which are derived towards a position of about 12" north- 
east from the target position of Bergin et al. (1999), are about 
3.8 ± 0.8, 9.3 ± 2.3, and 5.4 + 1.0 times higher, respectively. 
Moreover, the C ls O column density of ~ 5.8 + 0.2 x 10 15 cnr 2 
derived by Bergin et al. (1999) towards IRAS05399 implies the 
value N(C u O) ~ 1.6 ± 0.1 x 10 15 cnr 2 . The latter value is 
4.7 + 0.3 times higher than the value we have obtained. Besides 
a slightly different target position, these discrepancies are likely 
caused by the fact that Bergin et al. (1999) assumed optically 
thin emission, and the values n(H2) ~ 10 5 cnr 3 and 7\i n = 15 K 
in their non-LTE excitation analysis. 

The C n O column densities and abundances we have 
derived are comparable to those found by Bacmann et 
al. (2002) for their sample of seven prestellar cores, and 
by Schnee et al. (2007b) towards the starless core TMC- 
1C. The N2I-F column densities are also comparable to 
those determined towards several other low-mass starless 
cores and Class protostellar objects (Caselli et al. 2002b; 
I Crapsi et al.^0051 iRoberts & Millar 20071 iDanTel et al. 20071 
Emprechtinger et al. 2009| [Friesen et al. 2010a,b). On the other 
hand, the N 2 D + column densities we have obtained are generally 
larger than those derived for other low-mass cores in the studies 
mentioned above. Also, the H 13 CO + and DCO + column densi- 
ties and abundances we derive are higher than those obtained by 
Anderson et al. (1999) in the R CrA region, and by Frau et al. 
(2010) for a sample of starless cores in the Pipe Nebula. 

The HCO + abundance of ~ 1.4 x 10~ 8 we have derived to- 
wards IRAS05399 is relatively high. HCO + could, in principle, 
increase in abundance at later evolutionary stages, when the CO 
gas-phase abundance is enhanced, and the reaction + CO — > 
HCO + + H2 becomes efficient. However, the derived CO de- 
pletion factor in IRAS05399 is relatively large, fy ~ 3.2, and 
thus the high HCO + abundance may be the result of an alter- 
native production pathway, C + + H2O — > HCO + + H, viable 
in the shocks (Ra wlings et al. 2000[ IViti et al. 20021 . We note 
that IRAS05399 drives the HH92 jet, which lies close to the 
plane of the sky, and the presence of shock is therefore plausible 
(IGredel et al. 19921 |Bally et al. 2002| i. 



5.4. Depletion and deuteration 

The CO depletion factors we have derived are in the range /d ~ 
1.6 ± 0.5 - 10.8 + 2.2 (/ D = 3.6 + 1.2 for SMM4-LVC and 0.8 ± 
0.3 for Ori B9 N-LVC). Interestingly, the strongest depletion is 
observed towards the protostellar core SMM 3, and the lowest 
value is found in the starless core SMM 7. This could be caused 
by the fact that we do not observe exactly towards the submm 
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peak of SMM 7, whereas we probe the dense envelope of SMM 
3 (see Figs. [TJand[2]). The second lowest value, 1.9 + 0.3, 
is seen towards the starless core SMM 1. Bergin et al. (1999) 
mapped the IRAS05399/SMM 1-system in C ls O, CS, H 13 CO + , 
and DCO + , and found that the emission peaks coincide with the 
position of SMM 1, in agreement with our finding of low CO 
depletion. In general, the gas-phase CO abundance is expected 
to decrease during the prestellar phase of core evolution as a 
result of freeze out onto grain surfaces. During the protostellar 
phase, on the other hand, CO is expected to be released back into 
the gas phase via ice-mantle sublimation in the warmer envelope 
surrounding the protostar. 

Caselli et al. (2008) estimated the value / D (CO) = 3.6 to- 
wards a position at the edge of Ori B9 N using the data from 
Caselli & Myers (1995) and Harju et al. (2006). This is higher 
than the value 2.1 + 0.8 we have derived towards our target 
position near Ori B9 N (also at the core edge). For compari- 
son, Bacmann et al. (2002) studied the level of CO depletion 
in prestellar cores, and found the values in the range 4.5-15.5. 
These are comparable to the values we have obtained. Moreover, 
Bacmann et al. (2002) found a positive correlation between fu 
and n(H2) of the form /b oc n(H2) 4 ~ 8 . No correlation was 
found between and n(H2) in the present study. Emprechtinger 
et al. (2009) derived the values /b ~ 0.3 - 4.4 towards a sample 
of Class protostars. 

We note that in the depletion analysis we used the value 9.5 x 
10~ 5 for the undepleted CO abundance. However, this value is 
known to vary by a factor of ~ 2-3 between different star- 
forming regions. For example, Lacy et al. (1994) determined the 
CO abundance of ~ 2.7 x 10~ 4 towards NGC 2024 in Orion 
B. Adopting the latter value would result in ~ 2.8 times larger 
depletion factors. 

The N2I-F deuteration degree in the Orion B9 cores is found 
to be in the range tf D (N 2 H + ) =s 0.035 - 0.992. The extreme 
value of 0.992 measured towards SMM 1 is, to our knowledge, 
the highest level of deuteration in N2I-F observed so far. This 
suggests that the core is chemically highly evolved but is in 
contradiction with the low fy value which points towards a 
younger evolutionary stage. SMM 1 could have been affected 
by the outflow driven by IRAS05399, releasing CO into the gas 
phase and effectively resetting the chemical clock. In this case, 
the very high /?d(N2H + ) value would be remnant of the ear- 
lier CO-depleted phase, and not yet affected by the gas-phase 
CO which destroys N2FF and N2D + . An opposite situation was 
found by Crapsi et al. (2004) towards the chemically evolved 
dense core L1521F which harbours a very low luminosity object 
or VeLLO dBourke et al. 20061) . where f D ~ 15 but fl D (N 2 H + ) 
is only ~ 0. 1 . For comparison, Crapsi et al. (2005) derived the 
values of /?d(N2H + ) in the range < 0.02 - 0.44 for their sample 
of starless cores, and Daniel et al. (2007) derived comparable 
values of /fo(N2H + ) ~ 0.07 - 0.53 towards another sample of 
starless cores. Also, Roberts & Millar (2007) found fl D (N 2 H + ) 
values of < 0.01 - 0.31 for low-mass cores. Emprechtinger et 
al. (2009) found the values # D (N 2 H + ) < 0.029 - 0.271 towards 
Class objects. Pagani et al. (2009a) and Fontani et al. (2011) 
derived the value /?d(N2H + ) ~ 0.7 at the centre of the starless 
core LI 83 and towards the high-mass prestellar core candidate 
G034-G2 MM2, respectively, which were the highest fractiona- 
tions reported before the present work. 

The degree of deuterium fractionation in HCO + could be 
derived only towards IRAS05399 and SMM 4-LVC with the 
values fl D (HCO + ) 0.020 and 0.010, respectively. These 
are significantly lower than the corresponding /fo(N2H + ) val- 
ues, in agreement with the results by Roberts & Millar (2007) 



and Emprechtinger et al. (2009). Such a trend is believed to 
be caused by the role of CO, which is the parent species of 
HCO + and DCO + , but the main destroyer of H 2 D + , N 2 H + , 
and N2D + molecules. Therefore, molecular deuteration pro- 
ceeds most efficiently in regions where CO is depleted. In the 
warmer envelope layers where CO is not depleted, HCO + can 
have a relatively high abundance, resulting in a lower line-of- 
sight average value of .Rd(HCO + ) compared to that of N2H + 
(Emprechtinger et al. 20091. Note that we have employed the 
DCO + (4 - 3) transition, which is expected to trace the inner part 
of the high-density envelope, where heating by the embedded 
protostar may affect the deuterium chemistry [by lowering the 
7?d(HCO + ) ratio]. For comparison, J0rgensen et al. (2004) deter- 
mined similar 7?d(HCO + ) values of < 0.001 - 0.05 for a sample 
of 18 low-mass pre- and protostellar cores. 

Previous studies of low-mass starless and proto- 
stellar cores have found correlations between the de- 
gree of CO depletion and deuterium fractionation 
Bacmann et al. 20031 [J0rgensen et al. 2004[ [Crapsi et al. 2005 



|Emprechtinger et al. 2009 1. No correlation was found between 
fo and ^d(N2H + ) in the present study. However, the values of /d 
and fl D (N 2 H + ) for the protostellar cores SMM 3 and IRAS05399 
are in rough agreement with the finding of Emprechtinger et 
al. (2009), i.e., that deuteration in the envelopes of Class 
protostars increases as a function of fa. 

5.5. Fractional ionisation 

In Paper I, we determined the lower limits to x(e) of a few times 
10~ 8 , and upper limits of about six times 10~ 7 towards two target 
positions near the clump associated with IRAS05405. The lower 
limits to x(e) derived in the present work are comparable to those 
from Paper I. 

The standard relation between the electron abundance 
and the H2 number density is x(e) ~ 1.3 x 10~ 5 /i(H2)~ I/2 
(IMcKee 19 89). This is based on the pure cosmic -ray ionisa- 
tion with the rate 1.3 x 10 7 s _1 and includes no depletion of 
heavy elements. The standard relation yields the values x(e) ^ 
5.5 x 10~ 8 for IRAS05399, and 6.7 x 10" 8 for SMM 4-LVC. 
These are roughly comparable to the estimated lower limits to 
x(e). The values of £h 2 were found in Paper I to be ~ l-2xl0~ 16 
s . These resemble the value derived here towards SMM 4- 
LVC, but are an order of magnitude higher than obtained for 
IRAS05399. Instead, the value £h 2 ~ 2.6 x 10~ 17 s" 1 derived to- 
wards IRAS05399 is within a factor of two of the standard value. 
For comparision, observations of the Horsehead Nebula in Orion 
B by Goicoechea et al. (2009) could only be reproduced with 
<Th 2 =7.7±4.6x10~ 17 s- 1 . 

Caselli et al. (1998) determined x(e) in a sample of 24 low- 
mass cores consisting of both starless and protostellar objects. 
Their analysis was based on observations of CO, HCO + , and 
DCO + , and the resulting values were in the range 10~ 8 to 10~ 6 , 
bracketing the values for the Orion B9 cores. They argued that 
the variation in x(e) among the sources is due to variations in 
metal abundance and £h 2 \ the latter was found to span a range 
of two orders of magnitude between 10~ 18 - 10~ 16 s _1 . Some 
of this variation could be due to different cosmic -ray flux in the 
source regions. Williams et al. (1998) used observations of C I8 0, 
H 13 CO + , and DCO + to determine the values 10" 75 < x(e) < 
10" 6 5 in a similar sample of low-mass cores as Caselli et al. 
(1998), but using a slightly different analysis. They deduced a 
mean value of £h 2 = 5 x 10~ 17 s . Applying the same analysis as 
Williams et al. (1998), Bergin et al. (1999) found the ionisation 
levels of 10~ 7 ' 3 < x(e>) < 10~ 6y towards more massive cores 
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(in Orion) than to those studied by Williams et al.. Bergin et al. 
(1999) determined the fractional ionisation towards IRAS05399 
to lie in the range x(e) ~ 9.3 x 1(T 8 - 1.8 x 1(T 7 . The lower limit 
we have derived, x(e) ~ 1 .5 x 10~ 8 , is about six times lower than 
the corresponding value derived by Bergin et al. (1999). 

Recently, Padovani & Galli (2011) suggested that the mag- 
netic field threading the core affects the penetration of cos- 
mic rays, and can decrease the ionisation rate by a factor of 
~ 3 — 4. The values of (u 2 determined through observations (so 
far) would then underestimate the inter-core values by the above 
factor, making the ^ values more comparable with those mea- 
sured in diffuse clouds. 

5.6. Deuterated formaldehyde 

We have detected the ortho form of doubly deuterated formalde- 
hyde, D2CO, towards the prestellar core SMM 1 and the proto- 
stellar core IRAS05399. We derived the column densities and 
abundances of this molecule to be 1 - 2 x 10 12 cm~ 2 and ~ 
5 - 10 x 10~ n , respectively. We note that the high-temperature 
statistical ortho/para ratio of D2CO is 2, and appears to be simi- 
lar even in cold star-forming cores (see Roberts & Millar 2007). 

Using the IRAM 30-m telescope, Ceccarelli et al. (1998) 
detected D 2 CO towards IRAS 16293-2422. It was the first re- 
ported detection of D2CO towards a low-mass star-forming core. 
A lower limit to the D2CO column density they obtained, ~ 10 14 
cm -2 , is much higher than in our sources. Ceccarelli et al. (1998) 
speculated that such a high D2CO column density requires eva- 
poration of D2CO from the grain mantles, where it is expected 
to be formed during the prestellar phase. Later, Ceccarelli et 
al. (2001) demonstrated that heating by the central protostar in 
IRAS 16293-2422 is responsible for the mantle evaporation, and 
injection of D2CO into the gas phase. Loinard et al. (2002) de- 
tected D2CO towards 19 low-mass protostellar cores, and found 
the D 2 CO/H 2 CO abundance ratios of ~ 0.02 - 0.4. Bacmann 
et al. (2003) and Roberts & Millar (2007) found D 2 CO column 
densities of ~ 0.5 - 2.7 x 10 12 crrT 2 towards a sample of low- 
mass prestellar and protostellar cores. These column densities 
are similar to those we have derived. However, the D2CO abun- 
dances derived by Bacmann et al. (2003) are somewhat lower 
than those we have derived. Recently, Bergman et al. (2011) 
found the D2CO column densities of ~ 2 x 10 12 cirT 2 towards a 
few positions in p Oph A, which are very similar to our values; 
towards the D-peak of p Oph A, however, they derived a high 
value of ~ 3.2 x 10 13 cirT 2 . 

D2CO is expected to evaporate from the CO-rich grain 
mantles when the dust temperature exceeds about 25 K 
dCeccarelli et al. 200 II) . This would agree with the low CO 
depletion factor of 1.9 derived towards SMM 1, because 
the CO sublimation temperature is about 20 K (e.g., 
lAikawa et al. 20081 . It is uncertain, however, why we see higher 
degree of CO depletion towards the nearby protostellar core 
IRAS05399. Moreover, the gas kinetic temperature is only 11.9 
K in SMM 1 and 13.5 K in IRAS05399. It is possible, that 
the presence of gas-phase D2CO in these two sources is due to 
a non-thermal desorption mechanism (Roberts & Mi llar 2007b . 
Bergman et al. (2011) suggested that cosmic-ray heating and 
the formation energy of the newly formed species could play 
an important role in the release of D2CO from the grain mantles 
in starless cores. In addition to these two mechanisms, shocks 
caused by the protostellar outflow driven by IRAS05399 could 
be responsible in releasing D2CO from the icy grain mantles in 
the IRAS05399/SMM 1 -region. Towards IRAS05399 the emis- 
sion is perhaps dominated by the cool envelope whereas in SMM 



1 we might see a component which is interacting with the out- 
flow (causing the high gas-phase CO abundance). 

5.7. Fragmentation in the Orion B9 cores 

5.7.1 . SMM 6 - a fragmented prestellar core 

Figure [9] shows the SABOCA 350-//m image of the prestellar 
core SMM 6 overlaid with the LABOCA 870-^m contours. The 
core is filamentary in shape, and it is resolved into three to 
four subcondensations at 350 //m. The projected linear extent 
of core's long axis is about 1'9 or 0.25 pc, and the core's mass- 
per-length is Mn ne ^ 33 M0 pc~'. The length-to-width ratio in- 
creases from about 1 .9 at the NW end to about 7.4 at the SE end. 
The projected separation between the condensations is ~ 29" or 
~ 0.06 pc, where the sources 6c and 6d are treated as a single 
fragment. The measured separations should be taken as lower 
limits because of the possible projection effects. 

To analyse the fragmentation of SMM 6 in more detail, we 
calculate its thermal Jeans length from Aj = c 2 /(G£oX where 
c s = -sJIcftTkin/ /Amu is the isothermal sound speed (at 7\i n = 11 
K), G is the gravitational constant, and So = fimuN(H2) is the 
surface density; /j is the mean molecular weight per free particle 
(2.33 for He/H= 0.1), and N(H2) refers to the central column 
density for which we use the value ~ 10 22 cirT 2 [see Col. (3) 
of Table [§). The resulting Jeans length is ~ 0.05 pc, very si- 
milar to the observed separation between the condensations. In 
Paper II, we found from NH3 measurements that turbulent mo- 
tions within SMM 6 are subsonic, so their contribution to the 
effective sound speed would not increase the calculated Jeans 
length much. Assuming spherical geometry, the local Jeans mass 
of SMM 6 is M s = 4n/3 x <p)(/lj/2) 3 ~ 2.2 M . The cor- 
responding Jeans number is rtj = M/Mj ~ 4, which is similar to 
the number of observed subfragments. If the core substructure 
is caused by gravitational fragmentation, rtj is expected to give 
an approximate number of subfragments within the core (e.g., 
IRathborne et al. 20071 

The measured radii of the individual condensations, ~ 0.01 - 
0.02 pc, and their masses, ~ 0.1 - 0.2 Mq, are smaller than the 

Jeans lengths of A3 = yjncllGp ~ 0.05 - 0.06 pcQ and Jeans 
masses of ~ 0.8 - 0.9 M0 [assuming 7\; n = 1 1 K and densities 
from Col. (4) of Tabled . Because the subcondensations are pre- 
sumably colder than 1 1 K, their masses and densities are likely to 
be higher. For instance, using the dust temperature Tdust = 7.9 K 
derived earlier for SMM 6 would result in about 4.4 times higher 
masses and densities. Lower temperature would also cause the 
Jeans length and mass to be smaller. At 7.9 K, these would be 
~ 0.02 pc and ~ 0.2 - 0.3 Mq, respectively. Therefore, the ob- 
served condensation properties can well be comparable to the 
corresponding Jeans values. 

The above analysis suggests that thermal Jeans instability 
is the dominant process responsible for the core fragmentation. 
Moreover, the process of Jeans-fragmentation appears to have 
reached its final state at scale of the observed subcondensations. 
Thus, the condensations are potential sites to give birth to indi- 
vidual stars or small stellar systems. The condensations can grow 
in mass through (competitive) accretion from the parent core. In 
general, core fragmentation is believed to be the principal me- 
chanism for the formation of binary and multiple stellar systems 
(e.g.. lTohline^0^lGood"win et al. 20071 . 



12 For the condensations we utilise the /lj-formula which assumes 
spherical symmetry. For the filamentary parent core, we calculated the 
central Jeans length A] oc 1 /E - 
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Besides the Jeans analysis, it is interesting to examine 
whether SMM 6, owing to its filamentary shape, is unstable to 
axisymmetric perturbations. For an unmagnetised isothermal fi- 
lament the instability is reached if its Mi; ne value exceeds the 
critical equilibrium value of AfSJ 1 = 1c\jG (e.g. JOstriker 19641 
Inutsuka & Miyama 1992| >. For SMM 6, M™' * 18 M pc" 1 , 
and Mii ne /M™^ « 1.8. Thus, SMM 6 appears to be a ther- 
mally supercritical filament susceptible to fragmentation, in 
agreement with the detected substructure. We note that in the 
case of a magnetised molecular-cloud filament, MjP' differs by 
only a factor of order unity from that of an unmagnetised case 
(Fieg e&Pudritz 2000 1. 

The three-dimensional velocity dispersion in SMM 6 is 
<t 3d = 0.384 km s -1 (from the NH3 data in Paper II). This can be 
used to estimate the perturbation timescale, or the signal cross- 
ing time, T cross = D/(T3 D , where the diameter is D = 0.25 pc for 
SMM 6. The value of r cross is about 6.3 x 10 5 yr or ~ 2.8 times 
the free-fall timescale, Tff. This timescale is comparable to the 
lifetime of the prestellar phase of core evolution of a few times 
Tff (see Papers I and II and references therein). 

The observed substructure within SMM 6 shows that core 
fragmentation can take place already during the early prestellar 
phase of evolution, i.e., before the formation of embedded pro- 
tostar(s). Recently, Chen & Arce (2010) suggested that the three 
subcondensations in the prestellar core R CrA SMM 1A were 
formed through the fragmentation of the elongated parent core 
in the isothermal phase. It seems likely that SMM 6 has under- 
gone a similar fragmentation process. The universality of frag- 
mentation of the starless cores is unclear, however. For example, 
Schnee et al. (2010b) found that none of their 11 starless cores 
in Perseus are fragmented into smaller subunits. This is impor- 
tant knowledge when comparing the core mass function (CMF) 
to the stellar initial mass function (IMF), because the presence 
of substructure is believed to be one reason for the mass shift 
between the CMF and the IMF. 




Fig. 9. Detailed 350-yum SABOCA image of the fragmented 
prestellar core SMM 6 with power-law scaling. The contours are 
as in Fig. [2] The green plus sign indicates the position of our 
molecular-line observations. The colour-bar scale corresponds 
to Jy beam -1 . 



from the protostar position. Next to that, at the borderline of the 
3.3cr LABOCA 870-yum contour, there is another subcondensa- 
tion, SMM 3c. The projected separation between SMM 3 and 3b 
is about 0.08 pc. This is very close to the thermal Jeans length 
for the whole SMM 3 core of 0.07 pc. In SMM 7, the subconden- 
sation SMM 7b is at ~ 26" or 0.06 pc from the 'main' submm 
peak. Again, this is quite close to the local Jeans length of 0.09 
pc. It thus seems possible, that the substructure within SMM 3 
and 7 is caused by thermal Jeans fragmentation. We note that 
the SMM 3/3b-system is qualitatively similar to L1448 IRS2/2E, 
where there is a dust condensation (IRS 2E) next to a Class 
protostar; L1448 IRS2E is the most promising candidate of the 
"first hydrostatic core" detected so far dChen et al. 201 01) . 

SMM 4 is resolved into two fragments at 350 pm. The east- 
ern one, SMM 4b, is associated with a Spitzer source at 24 and 
70 //m. In Papers I and II, we proposed that N 2 H + could be 
depleted in the dense envelope of SMM 4. This was based on 
two observational results: i) the N2H + (1 - 0) map of Caselli & 
Myers (1994) shows no emission peaks near SMM 4; and ii) 
the N 2 H + (3 - 2) line was not detected towards SMM 4 at the 
systemic velocity ~ 9 km s _1 (see also Fig. [3]). The NH3 lines, 
on the other hand, can be seen at ~ 9 km s with an addi- 
tional velocity component at about 1.6 km s (Paper II). The 
present molecular-line observations, however, revealed a some- 
what surprising result: all the observed lines are at an LSR ve- 
locity of about 1.5-1.7 km s _1 , not at ~ 9 km s _1 . The absence 
of molecular- line emission at ~ 9 km s _1 is probably not due to 
chemistry, because it is not reasonable to think that all the ob- 
served species, in particular DCO + and N2D + , would have been 
depleted (e.g JLee et al. 20031 1. Instead, it seems that SMM 4 is 
a member of the low-velocity part of Orion B discussed in Paper 
II. Such a chance line-of-sight alignment is surprising, because 
the nearby sources SMM 5, IRAS05405, and Ori B9 N show 
line emission at ~ 9 km s _1 . The reason why we detected the 9- 
km s NH3 lines towards SMM 4 could be due to the large beam 
size (40") of the observations; NH3 emission could have been 
captured from the nearby cores by the beam. Also, the morpho- 
logy of the N2H + integrated intensity map of Caselli & Myers 
(1994) can be understood if the velocity range used to construct 
the map only covers emission near 9 km s _1 ; SMM 4 emits at 
lower velocity and does not show up on the map. The dust pro- 
perties of SMM 4, such as its SED, should be taken with caution 
as the core consists of subcondensations. It also seems possible 
that the emission from the clump near SMM 4 is 'contaminated' 
by an unrelated object(s) along the line of sight. 

6. Summary and conclusions 

We have carried out a (sub)mm study of dense cores in Orion B9. 
We used APEX to map the region at 350 fim, and to observe the 
transitions of C n O, H 13 CO + , DCO + , N 2 H\ and N 2 D + . These 
data were compared with our previous LABOCA 870-yum data. 
The principal aim of this study was to investigate the dust emis- 
sion of the cores near the peak of their SEDs, and the degrees 
of CO depletion, deuterium fractionation, and ionisation in the 
sources. Our primary results are summarised as follows: 



5.7.2. SMM 3, 4, and 7 

The SABOCA 350-//m images of SMM 3, 4, and 7 overlaid with 
870-/^m contours are shown in more detail in Fig. [TOj In SMM 
3, there is a subcondensation, we called SMM 3b, at about 36" 



1 . All the 870-//m cores within the boundaries of our SABOCA 
map were detected at 350 fim. The strongest 350-//m source 
in the region is SMM 3, a candidate Class-0 protostellar core. 

2. Four of the 870-jum cores, namely SMM 3, 4, 6, and 7, were 
resolved into at least two condensations at 350 fim. In partic- 
ular, the elongated prestellar core SMM 6 was resolved into 
three to four very low-mass subcondensations, showing that 
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Right ascension (J2000) 



Fig. 10. Zoom-in views of Fig. [2] with linear scaling towards 
SMM 3, 4, and 7. The box symbols mark the positions of the 
Spitzer 24-/mi point sources. The green plus signs show the po- 
sitions of our molecular-line observations. In each panel, the 
colour-bar scale corresponds to Jy beam -1 . 



core fragmentation can take place during the prestellar phase 
of evolution. In all cases, the origin of substructure can be 
explained by thermal Jeans-type fragmentation. 
3. The dust temperatures derived from the 350-to-870-/mi flux 
density ratio are very low, only 7d ust * 7.9*^2 - 10. 
K. The corresponding gas temperatures are typically a few 
kelvins higher. We also derived the submm dust emissivity 
spectral indices using the assumption T^st = T ga% , and they 
are in the range j3 « 0.5 ± 0.8 -1.8 + 0.6. The uncertainties in 
f5 are large, and within the errors the values are comparable 
to the fiducial value [3 — 2. 



We refined some of the protostellar SEDs presented in Paper 
I by adding the observed 350-/vm flux densities. No radical 
changes were found, and the sources IRAS 05405-01 17 and 
SMM 3 and 4 are classified as Class objects, in agreement 
with our previous results. 

The CO depletion factors were found to lie in the range 
fu ~ 1.6 ±0.5- 10.8+2.2. We found no systematic difference 
in /b between the starless and protostellar cores. In accor- 
dance with previous observations and theoretical predictions 
the most severe CO depletion is seen towards the core with 
highest average density, SMM 3. The degree of deuteration 
in N 2 H + was found to be in the range Af(N 2 D + )/./V(N 2 H + ) 
0.035 ± 0.006 - 0.992 ± 0.267. The N(DCO + )/N(HCO + ) ra- 
tio was found to be about 1-2%, comparable to those seen in 
other low-mass star-forming regions. 

The fractional ionisation could only be derived towards 
IRAS 05399-0121 and SMM 4-LVC with the lower limits of 
x(e) > 1.5 x 10~ 8 and > 6.1 x 10~ 8 , respectively. These val- 
ues are comparable to the fractional ionisations we derived 
earlier towards a few target positions near IRAS 05405-01 17 
and SMM 4. The cosmic-ray ionisation rate of H 2 implied 
by the derived lower x(e) limit is fn 2 ~ 2.6 x 10~ 17 s _1 
towards IRAS 05399-0121, and ~ 4.8" x 10~ 16 s" 1 towards 
SMM 4-LVC. The former value, which does not suffer from 
line-of-sight contamination, is within a factor of two of the 
'standard' value 1.3 x 10~ 17 s _1 . 

The highest degree of deuteration in N 2 H + , ~ 0.99, was de- 
rived towards the prestellar core SMM 1 . To our knowledge, 
this is the most extreme level of N 2 H + deuteration reported 
so far. We also detected D 2 CO emission towards SMM 1 
with an abundance of several times 10". Because D 2 CO 
is expected to be formed through the grain-surface chem- 
istry, its presence in the gas phase in SMM 1 could be due 
to shocks driven by the jet from IRAS 05399-0121, resulting 
in the release of D 2 CO from the grain mantles. This con- 
forms to the low CO depletion factor of 1.9 derived towards 
SMM 1 . The very high N 2 H + deuteration could be remnant 
of the earlier CO-depleted phase, and not yet affected by the 
destroying effect of gas-phase CO. 

It seems likely that the elongated clump associated with 
IRAS 05405-0117 and SMM 4 consists of physically inde- 
pendent objects that form a single clump only in projection 
along the line of sight. 



In Papers I and II we discussed the origin of dense cores 
in Orion B9. The region could have been affected by the feed- 
back from the nearby Ori OB lb subgroup of the Orion OBI as- 
sociation [resembling the process termed "cloud shuffling" by 
Elmegreen (1979)]. This feedback process may have played a 
role in sweeping up the gas into dense shells and filaments out 
of which the dense cores were later fragmented. This picture 
is supported by the kinematics of Orion B9 (low-velocity line- 
of-sight components). The present observational results suggest 
that the further fragmentation of cores into smaller condensa- 
tions is caused by the thermal Jeans instability. The observa- 
tions presented in this paper also reveal the intricate chemistry 
of the cores, such as heterogeneous depletion and deuteration 
among the cores. This suggests that the individual Orion-B9 
cores, though possibly formed collectively in the parent fila- 
ments, are evolving chemically at different rates. 
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